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CHAPTER I 


INTRODUCTION 


Recently, the indoor compact range has received much attention and 
is increasng in popularity as it rivals the outdoor range for antenna 
and scattering measurements. As the compact range performance improves, 
its use will continue to grow. An integral part of this system is the 
means of providing a uniform plane wave. This study presents a 
Cassegrain antenna feed system as a means to achieve a more uniform 
plane wave. 

Normally, a single parabolic reflector is used to generate this 
plane wave. Edge diffractions are the major drawback to this reflector 
system because they generate ripple on the desired uniform plane wave. 
One method used to reduce this ripple is through the use of serrated 
edges. The edges may also be rolled to reduce the ripple and provide 
structural strength as well. Large circular rolled edges provide 
greater ripple reduction but require additional structural support and 
are more costly. Elliptic edges are also used to control ripple and 
these are more effective than simple circular edges since there is more 
control over the shape such that a smaller ellipse can work as well as a 
larger circular edge. At the bottom of a parabolic dish section, 
serrated absorber patterns are often used to break up the field. Proper 
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tapering of the feed horn field pattern will also minimize the 
diffracted fields. 

Another possible alternative to the reflector system is the use of 
a lens antenna. But lens antennas are not being widely used due to 
several disadvantages. Although lens antennas are frequency dependent, 
tfie major disadvantage is in the construction of the antenna itself. 
Reflector antennas are much easier to design since only one surface 
needs to be considered. If made of natural dielectrics, lens antennas 
can be heavy and bulky, especially at lower frequencies. The 
homogeneity of the dielectric is also often in question. Besides the 
structural problems, lenses are also inherently lossy and reflections 
occur at both interfaces [1], Therefore, reflector type antennas are 
usually favored over lens antennas. 

The alternative considered in this study is the Cassegrain 
reflector system. Theoretically, the Cassegrain system offers better 
performance than a simple reflector system. The longer pathlengths in 
the Cassegrain system lead to a more uniform field in the plane of 
interest. The convenient location of the feed and supporting hardware 
is another advantage of this system. 

Initially, several disadvantages to the Cassegrain system are 
apparent. The addition of the subreflector increases system complexity 
both in terms of construction and performance analysis. The 
subreflector also gives rise to aperture blockage. The orientation of 
the feed now leads to spillover illuminating the target area as well as 
the rest of the room. Finally, the addition of the subreflector leads 
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to interactions between the two reflectors resulting in undesired field 
variations in the target area. These problems are considered as the 
Cassegrain system is designed and analyzed. 

The major design consideration in implementing the Cassegrain 
system is through the blending of the edges to improve performance as 
opposed to simply attaching circular or elliptical edges. The blending 
technique is a better method of providing the transition from one curve 
to another. Blending also reduces the junction diffracted field and 
therefore enhances performance. The tapering of the field is also 
controlled through the blending process. In fact, the blending concept 
is what allows the Cassegrain system to work as an effective source of a 
uniform plane wave. 
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CHAPTER II 


THEORY 


A. CASSEGRAIN REFLECTOR SYSTEM GEOMETRY 


The Cassegrain antenna system consists of a main reflector, 
subreflector, and feed. The main reflector is a parabolic curve while 
the subreflector has a hyperbolic contour. Two foci are present in 
this system: the real focal point located at the feed and the virtual 

focal point located at the focus of the parabola. To generate this 
system, two variables for each reflector must be specified. Seven 
variables are used to describe this system and three equations used to 
solve for the three remaining unknowns (Figure 2.1) are as follows: 


tan JUv = ±1 Dm 
7 T Fm 


( 2 . 1 ) 


l + l = 2 Fc 

tan <()V tan<j>r Us” , and (2.2) 


1 - sin( <j>v-(j)r)/2 - 2 Lv 

sin( 4>v+«t>r )/2 Fc . (2.3) 

Equation (2.1) applies to the main reflector while Equations (2.2) and 
(2.3) apply to the subreflector. The negative sign applies to the 
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Gregorian forms of the system. 


The parabolic main reflector is generated by 
Xm * 

4Fm . 

The hyperbolic subreflector is generated by 



where 

e = sin( <f>v+<f>r)/2 
' sT’n( <}>v-<|>r)/2 


a 


Fc 

Te , and 


(2.4) 


(2.5) 


( 2 . 6 ) 

(2.7) 


b = a/e2-i 


( 2 . 8 ) 


These equations govern the classical Cassegrain s>y:>teiiu Using thei>e 
equations, many variations of the basic system may be formed. For this 
study, the basic system of Figure 2.1 is sufficient, though one 
variation is considered. The Gregorian form occurs when the focus of 
the main reflector moves between the two reflectors (see Figure 2.2). 

In this case the contour of the subreflector is elliptical. The 
negative sign must be used in equation (2.1) since <f>v is now negative. 
Otherwise, the equations remain the same. 

In later analysis, the use of the virtual feed is made. In this 
concept, the real feed and subreflector are replaced by a virtual feed 
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Figure 



.1 Geometry of Cassegrain system. 



Figure 2.2 Gregorian form. 


at the focal point of the main reflector. The system is now a single 
reflector design, and this concept is useful when designing and 
analyzing the main reflector. 

Another useful concept is that of the equivalent parabola. As is 
seen in Figure 2.3, the Cassegrain system is replaced by an equivalent 
surface which has a parabolic contour as demonstrated by Hannan [2]. 
Using this concept, the feed remains unchanged, and ray optics are used 
to determine the surface as the locus of incoming rays intersecting the 
rays converging to the real feed. Then the Cassegrain system has been 
replaced by an equivalent single reflector system. The following 
equations show the relationship between the two systems: 


= tan Ur 
4 te 7 

(2.9) 

Xe = Ye 2 


4Fe . and 

(2.10) 

+ Fe = tan<j>v/2 = L r _ e +l 


Fm tan <t>r/2 e-1 . 

(2.11) 


Again, the negative sign applies to the Gregorian forms. Equation (2.9) 
generates the equivalent focal length. Equation (2.10) generates the 
equivalent parabolic surface, and Equation (2.11) provides several 
expressions for comparing the focal lengths. For the classical 
Cassegrain system, Fe/Fm is greater than one. It is apparent that a 
Cassegrain system has a much smaller focal length but can be equivalent 
to a single reflector system of much larger focal length. It is this 
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idea that favors the Cassegrain system over the single reflector system. 
When working in a restricted area, such as a compact range, the shorter 
Cassegrain system is favored over the longer equivalent single reflector 
system [2], 

B. MOMENT METHOD ANALYSIS 

Several analysis techniques are used when studying the Cassegrain 
system. The first that will be described is the moment method theory. 
Only the two-dimensional part of this theory is outlined [3], By using 
the reaction concept of Rumsey [4], a moment method solution may be 
obtained. In Figure (2.4) a scattering problem is presented. The 
source electrical and magnetic currents (J-j ,M-j ) generate electric and 
magnetic fields (E,H) in the presence of the scatterer which is a 
conducting body in free space. 

The surface-equivalence theorem of Schellkunoff [5] gives an 
equivalent problem by replacing the scatterer by the following surface 
current densities: 

J s = n x H , and (2.12) 

M s = E x n (2.13) 

a 

with n being the outward normal to the surface. It is self-evident that 
the source currents (J n * ,M-j ) generate the incident fields (E-j ,Hi ) in the 
free space. The scattered fields are given by 

E s = E - Ei, and (2.14) 
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-► -► + «► 

(a) The source (Ji ,M-j ) generates the field (E,H) with 
the scatterer present. 



+ -► 

(b) The interior field vanishes when the currents (J S ,M S ) 
are introduced on the surface of the scatterer. 


/ 

/ 

ll- 

\ 




\ 

I * N 
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(c) The exterior scattered field may be generated by (J S ,M S ) 
in free space. 


Figure 2.4 Scattering problems 


H s = H - Hi . 


(2.15) 


The surface currents generate the scattered fields (E S ,H S ) outside the 
body and ( -E^ , -H-j ) inside the body. 

An electric test source J m is now placed within the region of the 
scatterer (see Figure 2.5). Because there is no field in this region, 
the reaction of this test source with the field produced by the other 
sources is zero. By reciprocity, this reaction must be equivalent to 
the reaction of the other sources with the field produced by the test 
source. Then, one finds that 

//( J s* E m- M s* H m) ds + J7J( J i * E m" M i * H m) dv = 0 • (2.16) 

This equation is the basis of this solution and approach is the 
"zero-reaction theorem" of Rumsey [4], 

Next, the surface current distributions (J S ,M S ) need to be 
determined. These currents are constructed of finite series with N 
unknown coefficients. For this problem, the scattering body is a 
perfect electrical conductor so M s vanishes. Since only the 
two-dimensional case is being considered, J s is a function of the 
position z around the contour of the body. Also consider a magnetic 
line source and TE Z polarization so that is zero. Then the integral 
equation reduces such that 

J J s‘ E m d* = // Mi -H m ds . (2.17) 

c 
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lit J8.»\ 



\ FREE JFACE. 

s. y 



FREE SPACE 


Figure 2.5 Placement of test source. 

The electric currents may now be represented by 
N 

J S U) = I W*) (2.18) 

n=l 

where I n are complex constants, and samples of J s (t) and J n (£) are the 
basis functions. The basis functions as well as the test source have 
unit current density at their terminals. Substituting this series into 
the integral equation yields 
N 

I In^mn = ^m with m - 1,2,3,...N (2.19) 

n=l 
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where 


^mn = " In = "Jm'W*') *^nd4» ar *d (2.20) 

s ”//l^1 *^m^ s s Jm'W*) ‘^idt (2.21) 

with the integrations are done over the non-zero range of the 
integrands. 

When solving these expressions on a computer, it is advantageous 
for the impedance matrix, Z mn , to be symmetrical . Also, the test 
sources, J m , should be the same size, shape, and functional form as the 
basis functions J n . This will allow some of the integrals to be solved 
in closed form and yield readily solvable simultaneous linear equations. 
In addition, the test sources are placed a distance 5 from the surface 
where 6 tends toward zero in the limiting case of the integrals. In 
this case, sinusoidal strip dipoles are chosen as the basis functions 
(see Figure 2.6). This planar strip dipole extends infinitely in the 
z-di recti on and has a surface current density given by [9] 

j = x sin[k(x-xi)] 

sinLk(X 2 -xi)J (2.22) 

for xj < x < X 2 > and 

A 

j = x sin[k(x 3 -x)] 

sin|_k(x 3 -x 2 )J (2.23) 

for X 2 < x < x 3 . Similarly, for the strip v-dipole in Figure 2.7 the 

surface current density is given by 

j = _ s sin[k(ti-t)] 

sin(ksi) (2.24) 
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Figure 






(a) A planar strip dipole with edges at xj and X3 and 
terminals at xg. 


y 


L 


(b) 



The current -density distribution J on the sinusoidal 
strip dipole. 


2.6 Sinusoidal strip dipoles for basis functions. 
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on arm s, and 


j = t sin[k(tj-t)] 

sin(ktj) (2.25) 

on arm t. It is evident that s and t are perpendicular to the z-axis. 

In both cases, the current density goes to zero at the endpoints and is 
unity at the center terminals. Also, a slope discontinuity is present 
at the center terminals. Of course the v-dipole reduces to the planar 
case when is equal to 180°. 


The field of the strip dipole may be obtained from the 
superposition of two strip monopoles considered to have a common 
endpoint (see Figure 2.8). The field for the strip monopole is known 
from reference [6], Now the calculation may begin for elements of the 
impedance matrix as well as the V m elements of the excitation column. 

First consider an open or closed perfectly conducting polygon 
cylinder (see Figure 2.9). For this open cylinder, surface currents 
flow on both sides of the conductor, and the surface current density is 
given by J s . A magnetic line source M-j is present and Ij and I 2 are the 
current densities at the corners of the polygon. Two strip dipole mode 
currents may now be defined on the conductor. The first extends from 
point 0 to point 2 with terminals at point 1. The second extends from 
point 1 to point 3 with terminals at point 2. Each mode has a 
sinusoidal current distribution as described earlier. Now the current 
density J s is the superposition of these two modes with weightings of Ij 
and I 2 . Then J s is a piecewise sinusoidal expansion with unknown 
constants, Ii and I 2 . Since the polygon is a perfect conductor, the 
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3 


• M,‘ 



(a) Perfectly conducting polygon cylinder with parallel 
magnetic line source M-j . 



(b) Electric test probes 1 and 2 are moved to the conducting 
surface. 

Figure 2.9 Probing of a perfectly conducting polygon cylinder. 
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tangential electric field is zero on the surface. So, if an electric 
test probe is moved along the conducting surface, the open circuit 
voltage at its terminals can be examined. For N different current 
samples, N probing tests are done. The probes may be real (thin wire 
v-dipoles) or hypothetical (electric line sources or strip dipoles). 

Then the currents I n are adjusted until all the probes read zero. 
Finally, as N increases this stationary solution for the currents 
approaches the rigorous solution. The mutual impedance between the m^h 
test probe and the n^* 1 current mode is Z^. The probe sums all the 
voltage contributions from J s and M-j and this result must be zero, 
resulting once again in Equations (2.19) to (2.21). 

Finally, the simultaneous linear equations are solved using linear 
algebra techniques, and the current distribution J s is known. The 
scattered fields (E S ,H S ) may then be found. J.H. Richmond [6] provided 
this theory, method, and appropriate computer programs. Using duality, 
the TM polarization case could also be solved. 

C. UNIFORM GEOMETRICAL THEORY OF DIFFRACTION ANALYSIS 

Another analysis technique used is that of the Uniform Geometrical 
Theory of Diffraction. Again, only the two-dimensional part of this 
theory is outlined. 

For the Cassegrain system three basic field components are examined 
(see Figure 2.10). These are the incident, reflected and diffracted 
fields. The total field is then given by 

IjTOTAL „ ijINC + (jREF + „DIF . . 
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where U represents an electric scalar field for the electric line source 
case and a magnetic scalar field for the magnetic line source case. The 
incident field is given by 


ijInc 



Regions I and II, and 


Region III 


(2.27) 


whereas is the distance between the source and receiver. Note that C 
is a complex constant. The reflected field is given by 


.REF . 


±C 


^-jkpr 

~7oT~ 


Region I, and 
Regions II and III 


(2.28) 


where p r is the distance from the image of the source to the receiver 
and the positive sign is used for the magnetic line source case and the 
minus sign for the electric line source case. To simplify calculations 
the magnetic line source is used when analyzing the Cassegrain system. 
The diffracted field is given by 


u DIFF = 

D 

p' p , 4>— <f> ' ,n 

± D 

p‘ p » »n 



_p'+p _ 


_p‘+p _ 


-jkp' -jkp 
e 

/p 7 " /p* 

(2.29) 


with the plus and minus signs for the magnetic and electric line source 
cases, respectively. The term with is associated with the incident 
shadow boundary while the term with <}>+<}>' is associated with the 
reflected shadow boundary. 
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Now consider a curved conducting strip as will be found in the 
Cassegrain system (see Figure 2.11). The incident field does not change 
but the reflected field is now given by 


U REF * 



/ 


pc (Or) 
PcTQrT+Pt 


-jkpi -jkp r 
C l e 

/p? 


Region I, and 


(2.30) 


Regions II and III 


with the calculation of p c (Qr) needed. This caustic distance varies 
with the reflection point Qr and is given by 

1 =1 + 2 

Pc(Qr) Pi Rc cos e i (2.31) 


where R c is the radius of curvature of the surface at the reflection 
point, Qr. The diffracted field is given by 


u DIFF = 

D 

p' p , 4>- 4»' ,n 

± D 

Pc'p .(fr+^.n 



_p'+p 


_pc'+p _ 


-jkp' -jkp 
q e_ e 

'/p r /p* 


(2.32) 


I 

where p c is the caustic distance p c (Qe) associated with the reflected 
ray from the edge and n=2 such that 


-e“ jir/4 F(KLa( $)) 
D(L,B,n=2) = - 2/2 ^ cos (3/2) 


(2.33) 
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00 


F(KLa( $)) = 2j/KCaT0Te dkLa(P) 


/ 


✓KLa'(e) 


“J t2 

e dx. 


and 


(2.34) 


a(3) = 2 cos 2 (3/2) . 


(2.35) 


The diffracted field given is sufficient for the knife edge case 
but the general form of the diffraction coefficient is also needed. So 
more generally, the diffracted field is expressed by 


U°IFF = d s U 1 (Q e ) 
H 



(2.36) 


where U^(Q E ) is the incident field on the edge, and 

D s U»*'.n) = 2n/2irk x 
H 

[cot ( *+(♦-♦*) ) F [KL i a + ( <j>* )] + cot (*-(♦-♦') ) F [Kl^a’f 4»-4>* ) D+ 


(cot (llUltH) F [KL rn a + ( 4>+4>' ) ] + cot ("-(♦+♦*) ) FCKl/VU+V )]}] 
Zn Zn 


(2.37) 


where 



a ± (B) = 2cos 2 (2nirN ± -B) 


(2.39) 


2 

Note that N* are Integers most nearly satisfying the following 
expressions: 

2irnN + -(B) = it, and (2.40) 

2imN"-(3) = -ir. (2.41) 

These expressions are valid for both the soft and hard diffraction 
coefficients but only the hard case is used here with a magnetic line 
source. The transition function F(x) involves a Fresnel integral, and a 
plot of F(x) is shown in Figure 2.12. 

The diffraction coefficient may also be written as 

D s = D(Ll,Ll,M',n) + D(L rn ,L r0 , ,n) (2.42) 

H 


where 

D(Li,l_ 2 ,<t>±<|>',n) = [cotf TT+Utifr 1 ) ) p ( KLia + ( ) ) + 

-j ir/ 4 

cot f ir-( (frtift 1 ) F(KL2a-(<J>±<t>'))] 

2nvCT • (2.43) 

Figure 2.13 shows a more generalized structure. Depending on the 
positioning of the line source, reflected fields may emanate from both 
surfaces. So two reflection shadow boundaries may exist, and hence the 
reason D(L rn ,L ro , ,n ) is composed of two terms. The first term is 
associated with the n face boundary and the second with the o face 
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Figure 2.13 One face of a general wedge structure is illuminated. 
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boundary. The o face is defined as the face where the angles <$> and <p' 
are measured. The n face is the opposing surface. Also, the range 
parameters are given by 

L rn = PcP , and (2.44) 

pjpp 


o 


l_ro _ p c p 


it— 

p c +p 


(2.45) 


u n 0 

where p c and p c are the caustic distances for the reflected waves 
emanating from the edge for the n and o faces, respectively. Similarly, 
two incident shadow boundaries may exist, and two terms are also 
associated with these. In this case, the range parameter is given by 

L i = p'p 

7^ (2.46) 


This completes the basic theory and more detailed analysib may be found 
in the class notes by Burnside [7] for microwave optics. In later 
chapters some additional details of the theory are needed, and they are 
presented when needed. 

One final technique used in this analysis is physical optics. This 
involves rather simple analysis and will be described later when 
actually implemented. 

The three analysis techniques may now be compared to see the 
advantages and disadvantages in each for analyzing the Cassegrain 
system. The Uniform Geometrical Theory of Diffraction (UTD) requires 
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analysis of the geometry which may or may not be easily implemented. 
Including enough terms to accurately predict performance may be 
difficult but if possible, UTD provides results with very little 
computation time required. The UTD is also well suited for large 
electrical objects such as the Cassegrain system. If results are 
consistent with other techniques, UTD may be used as a valuable design 
tool because of its high frequency capability and ease of computation. 
The physical optics technique is also easy to implement and its results 
may be easily compared with the other techniques. Physical optics is an 
approximation though, and this is a limitation. Finally, the moment 
method technique provides the greatest accuracy but at the expense of 
ease of computation. Much computational work must be done so moment 
method results require much time and space on a computer. The moment 
method is also limited by object size. In the case of the large 
Cassegrain system, this limits the upper frequency that may be examined. 
Since the moment method provides accurate results, it may be compared 
with UTD to see what field components are dominant. So the faster UTD 
may be used to initialize a design and the slower moment method to 
finalize it. The moment method will also give more accurate results at 
lower frequencies; whereas, UTD may be used to predict the high 
frequency behavior. 
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CHAPTER III 


CASSEGRAIN SYSTEM CONSIDERATIONS 




A simple single parabolic reflector is considered initially. A 
line source is placed in the presence of half a parabola (see Figure 
3.1) whose focal length is specified by f. The reflected field and the 
edge diffracted field are examined along a line parallel to the x-axis 
and a distance f from the origin. All path lengths from the focal point 
to the reflector and onto this line are equidistant and have a value of 
2f implying constant phase across the plane for the reflected field. 

Now the reflected field in Figure 3.2 is given by the geometrical optics 
expression such that 


..RFF 

U”~ 


± 7 pc 

/ p c +p r 


-jkpi 

* n 

7 W 


-jkpr 

e 


13. ij 


The caustic distance is given by 


!_ = !_+ 2 

Pc Pi ft c cos(ir-9-j) (3.2) 


where 


R 


c 


1 

CURVATURE' 


and 


(3.3) 


CURVATURE = y " • 

(l+(y' ) 2 ) 3/2 


(3.4) 


29 




Figure 3.1 Parabolic reflector. 



Figure 3.2 Reflected field. 
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Solving these expressions for R c , one finds that 
R c = |2f(l+x 2 /(4f 2 )) 3/2 | . 

In addition, one obtains the following: 
cos ( 2 9-j ) = p c / pi 
cos 2 0-j = (1 + cos(20-j )/2 
cos 0^ = /[l+p r /pi )/2 . 

cos ( it— 0-j ) * -cos0i , and 

p r + pi = 2f . 

Substituting thse into Equation (3.2), one obtains 


(3.5) 


(3.6) 

(3.7) 

(3.8) 


Now 



1 

1 . 

f [l+x2 13/2ff jl /2 
k 4f2 J 


(3.9) 


Pr = f - x2/(4f) 
or 

pr= i . x 2 (3.10) 

f 4f 2 • 

Substituting this result into Equation (3.9), one finds as expected 
that 


p c = «° , and 


(3.11) 
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J REF = ± C i 


-jk(pi+p r ) 




vith 


p r = f - x^/(4f) , and 

Pi = A 2 + p 2 • 

Now U r ^F i s known as a function of x. 

The edge diffracted field in Figure 3.3 is given by 

p'p -jkp' 

U DIFF = [D(£!_£_, <j>- < j>* » n ) ± Dli^H'-njlc* 




Since n * 2 in this case, one finds that 

-j if/ 4 


D(L,B,n-2) = “ e 

F ( kLa ( e) ) 


2/2itT 

cos (0/2) 

jkLa(0) « -jx 2 

(3.16) 

F(kLa(0) = 2j/kLa( e) 

e / e dx, and 

/FTaTe)' 

(3.17) 

a(e) = 2cos 2 ( 0/2) . 


(3.18) 


From Figure 3.3, one obtains 



p 2 = x 2 +(p') 2 - 2xp'cos0" 

6" = ir/2 - 0' 

0' = sin -1 (a/p') 
cos ( ir/2 - 0') = sin 0' 
sin0' = a/p' , and 

p ■ (x 2 + ( p' ) 2 - 2xa)!/ 2 . (3.21) 

Now (f)' and <)> need to be determined. From Figure 3.4, one obtains the 
following: 


sin(20 1 ) = a/p' 

0 1 - 1/2 sin~l(a/p' ) (3.22) 

♦ ' * V2 - 0 1 (3.23) 
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x 2 = ( p ')2 + p 2 _ 2p' pcosy 


, . cos-^ -X^P') 2 U? ) ( 3.2 4) 

<(> = + y (3.25) 

0 " = + - ♦' 

0" = Y » and (3.26) 

0 + = 4> + <t>' . 

Substituting Equations (3.25) and (3.23) into (3.22), one finds that 

0 + = ir - sin"l(a/p') + y . (3.27) 



Figure 3.4 Edge diffacted field. 
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The range parameters are given by 


L" = p'p . and (3.28) 

p +P 

I 

L + = p cp = p (3.29) 

P +P 

c 

■ 

since p c = <®. 

Now the edge diffracted field is known from 


n -j*/4 -jk( p'+p) 

U DIFF = If ( F(kL-a(r)) ± F(kL+a(g+)) l c e . 

2/£irk‘ cos(g"/2) cos(g + /2) /pV* (3.30) 


Then the total field is given by U REF + U DIFF in Equations (3.1) and 
(3.30). A typical plot is given in Figure 3.5. 
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Now the addition of a rolled edge is made to the parabola to reduce 
the ripple generated by the diffracted field (see Figure 3.6). The 
diffracted field at the junction must now be recalculated. More 
generally, one obtains that 


-jkp 

U DIFF = D S uMQe) f 

H /T 

-jkp' 

U 1 (Qr) = C e 

11 ST 


(3.31) 

(3.32) 


D S = D(U,Li,4»-<|) , ,n) + D(Lrn, L ro,*f*', n ), and (3.33) 

H 

D(L 1 ,L 2 ,<|.±^',n) = [cot [ tt+( (t>t<j>* ) )F(kLia + ( )) 

-j w/4 

+ cot f tt-( ) )F(kL?a~( (frill) 1 ) ) 1 e . (3.34) 

^ 2n/2irk' 


The distance given by p and p 1 have been previously calculated but the 

two terms of D$ need to be considered. Looking at D(L"* .L 1 , 4>-<f>'»n), one 
H 

finds that 

. . -J*/4 

D ( L 1 jL 1 , ,n) = -e [cot [ lli±y± )V (kC a { ♦-»■)) 

2n /ZiT 2n 


+ cot f "-( ) )F(kL i a'f <p- <(,')) ] 

(3.35) 

n = 1 

(3.36) 

L 1 = p' p 

(3.37) 


p'+p 
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Figure 3.6 Rolled edge addition, 
and <f> and <t>' are also known. Recall that 
a ± ( 0) = 2cos 2 r 2mrN*-(3) ) 

where N* are integers which most nearly satisfy the following: 

2irnN + -($)*» , and 

2imN" - ( 3 ) = -ir . 

But n=l, so 

a ± (3) a 2 cos 2 (ttN ± - 8/2) , or 
a ± ( 8) = 2cos 2 (8/2) 

which implies that 
a + (8) = a"(3) . 


(3.38) 
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Inserting Equation (3.38) into (3.35) yields 
. . -J*/4 

D(L\L\<MM) = lf {cot (ir/2 +U-*')/2) 

2/7nF 

+ cot (it/ 2 -(♦-♦' )/2)} F(kL i a + ( ) . 
From the geometry, one finds that 


COt(ir/2 + a) + cot (n/2 - a) 


COS(tt/ 2 + a) + COS(tt/2 - 
sin(ir/2 + a) sin(ir/2 - 


= -sina + sina 
cosg cosg 


= 0 . 

Thus, the incident shadow boundary terms are given by 

D(L i jL 1 »1) = 0 

Now looking at the remaining terms, one obtains 

-j*/4 

D(L rn ,L ro , 4> +< i ) ' ,1) =— [cot (rr/2 + ( 4>+4> *)/2)F(kL rn a 

2/2^F 

+ cot(ir/2 - ( <J>+4> ' ) / 2 ) F ( kL ro a” ( * ) )]. 


(3 

«) 


(3 

*(♦+♦')) 

(3 


But as before. 


Now 


L ro = p?p_ = p (3.43) 

P c + p 

0 

since p c is infinity from previous calculations. 

For the n face, one obtains that 

L rn = Pc°P , and (3.44) 

n x 

Pc + P 

1 = 1 _ + 2 . 

p" p' RccosQ 1 (3.45) 

But p' is known and R c is given leaving 

e 1 = 1/2 sin” 1 (a/p l ) 
or 

cose 1 = cosCl/2 sin" 1 (a/p')] . (3.46) 

n 

Note that p c is known from Equation (3.45) and L rn from Equation (3.44). 
Finally 4>+ <f>' is given by Equation (3.27) and the diffraction coefficient 
is known from Equation (3.41). The diffracted field is then obtained 
using Equation (3.31). Again the reflected field and recalculated 
diffracted field are summed for various radii of curvature. In Figure 
3.7, the reduction in ripple is evident with increasing radii of 
curvature. The discontinuity at the edge is also apparent resulting 
from the absence of the edge reflected field. 
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Figure 3.7 (Continued). 


Now consider the Cassegrain system. The governing equations have 
previously been described in the theory section. Again, the reflected 
field is analyzed first (see Figure 3.8) and is given by 


jREF = 

1 / 


~ j k Pi • j k pi 

c e e 


cl^rl 


(3.47) 


The reflected field caustic position is calculated easily from the 
following geometrical considerations: 


Pci * ( Lv ‘ x S ) 2 + *s 2 (3.48) 
p i 2 = y s 2 + ( Fc ’Lv+x s ) 2 , and (3.49) 
p rl 2 = (V y s )2 + ( p r2 " ( Lv -* s ) )2 • (3.50) 
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Now it is necessary to relate y s in terms of y m , the desired field 
point. First, one obtains that 


2 

P r2 = Fm - y m /( 4Fm ) 

x s = a C /1 + (y s /b ) 2 - 1 ] , and 

tan a = Pr2 = *-v-x s 

y m ys 


Substituting Equations (3.51) and (3.52) into (3.53) yields 


'1 + 


ys - 


Lv + a - a' 


ier 


Lv 2 + 2aLv 

(hCj ) 2 


Figure 3.8 Reflected field 



So now P c j, Pi, and P r i are in terms of y m , and UjjEF may be 
calculated from Equation (3.47). Then, the reflected field along the 
observation line is given by 


U* EF = u FEF e ~^ kPr2 


(3.56) 


since Note that p r 2 is given by Equation (3.51). 

The diffracted field from the main reflector edge is computed now. 
As with the parabolic problem, one finds that 


U? IFF = U? FF L. 1 [Df -(^PC1)P ,<j»+i»',n) + P( Pc2P „fH>',n)] 
/r E . E PC2+P 

p rl +p cl +p 


'l 


’IE 


(3.57) 


where 



(3.58) 


E 

Pr- = 


C *y(Lv-Xs) 2 +(Ds/2)2 


= a bj\+ 


(^i) -1] 


s "\/ds/ 2) 2 + (Fc-Lv+Xs) 2 


(3.59) 


(3.60) 


(3.61) 
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(Dm/2-Ds/2) + (p^-(Lv-x‘)r , and 


p rl = ' 

Pp 2 = Fm - (Dm/2) 2 . 


Again, the diffraction coefficient consists of terms as follows 


D(L, B,n=2) 


-jtt/4 

F(kLa(g)) 

2 rnt cos (6/2) 


and 

a( 6) = 2cos 2 ( 6/2) . 

From parabolic edge calculations, one obtains that 


’ * C + ( (, rl + 4> 2 - *r» Dm > 1/2 

'tr*< 

6" ■ ♦ - 

6" = cos -1 [ -ym^Prl+Pcl) 2 + p 2 ] 
2(p rl +p cl )p 



L + = P 


(3.68) 


since P c 2 = «. So finally, the diffracted field is given by 


-jkp -jt/4 

U? IFF = U? FF 1 • if fF(kL-a(e-)) + F(kL+a(0+))1 . (3.69) 

2/jT 2/2irk cos (0-/2) cos(0 + /2) 


The sum of U^F and ipIFF i s shown in Figure 3.9 for one example. 

2 1 

Another major field component is the diffracted reflected field 
from the subreflector edge (see Figure 3.10). This field is given by 


.D IFF 


■/ 


-jkp n ' -jkpj £ 

HO c ! ! [nf pj p l ,<fr-4>',2i 

Pc + P2 /T /pT E* 

/pi K 1 p i +p l 


* DC£?L.*+r.2)]}e 


-jkp. 


O^+D- 

' C ’ i 


( ir\\ 

V ^ v / 


with 


D(L, 0,n=2) 


-jn/4 

-e F(kLa( 0) ) 

2/2ttF cos (0/2) 


and 


a(0) = 2cos 2 (0/2) . 


(3.16) 


(3.18) 
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The geometry for this field analysis is given in Figure 3.11. The 
needed relationships between the various parameters are given as 
f ol 1 ows : 


X m = (Dm/2) /(4Fm) 

(3.71) 

s 2 = (Xm - (Fm-Fc) 2 ) 2 + (Dm/2) 2 

(3.72) 

9 1 = 1/2 cos _1 [( p i ) 2+ ( p rl) 2 " s2 ] 

, E E 
2p i p rl 

(3.73) 

♦' 9 1 

(3.74) 

(s') 2 - y™ + (Em-Fc-y^/ (4Fm) 2 

(3.75) 

6- • *-*• > CO,' 1 <P1) 2 +(P1) 2 -(S‘) Z 

2p 1 p l 

(3.76) 

3 + = <t>+<t>* = 3~ +ir-20. 

/ o ~i -t \ 

) 

E 

L" = p i p l , and 
E 

p i +p l 

(3.78) 

. E 

L + = PcPl . 

p c +p l 

(3.79) 

p c in Figure 3.12 is given by 


-L = 1_ + 2 

Pc PI R c cos(Tr-9 i ) 

(3.80) 
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and 


2 2 3/2 

2Fm(l+y m /(4Fm ) 


R 



(3.81) 


from Equation (3.5). The reflection point is found from an iterative 

routine given in Appendix A. Now U^IFF -j s known from Equation (3.70). 

The sum of ipIFF and ipIFF -j s given in Figure 3.13 as a sample 

2 12 

geometry. 

The discontinuity in Figure 3.13 may be alleviated somewhat by 
extending the subreflector and reducing doubly diffracted fields (see 
Figure 3.14). The subreflector size is now determined by <|> s ; whereas, 
Dg may be found from 


TAN 4 * = De/ 2 
Lv-xs 


D E /2 


(3.82) 


Inserting the system parameters into this equation and solving for Df 
yields 


D F = (TAN <i> ) [a+Lv±/a 2 +(2aLv+Lv z )( (a z TAN z <t>s)/b z )] (3.83) 

t s 2 2 

1 _ a TAN 

2 2b 2 

with the minus, sign giving the desired solution. 

The diffracted reflected field may now be recalculated (see Figure 
3.15) using 


DEED 

x 

s 



(3.84) 


DEED 

P. 

1 


= / ( d E/2) 2 + (Fc-Lv+x s ° EED ) 2 


(3.85) 
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Figure 3.14 Subreflector extension 



Figure 3.15 Subreflector diffracted field. 


p° EED - /(D E /2) 2 + ( L v - x ^ EED ) 2 f 


and 


(3.86) 


DEED 

____ -jkp-i -jkpi DEED 

J° IFF ■ /3T( c * * [P( * 1,1 ,♦-♦■.2) 

* PC+P2 /DEED 


DEED 

+ D(J* 2) ]}e 

DEED , 
p c +p l 

In this case, one finds that 


1 

-jkp. 


DEED A 
p i +P 1 


(3.87) 


29^ * ir - a 


where 


a = cos 


-1 


. DEED 2 , DEED 2 2 

[( Pi ) +(p c ) -Fc j 

0 DEED DEED 

2p, o c 


(3.88) 
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In addition. 


♦' = w /2 - e. # 

Again, g“ is given by Equation (3.76); whereas. 


(3.90) 


B + = $+<{,' = g% a . 


(3.91) 


The sum of U^, U j* FF and is again shown in Figure 3.16. The 

reduction in the discontinuity is obvious. 

To reduce the ripple, a rolled edge is now attached to the main 

reflector (see Figure 3.17). The analysis is similar to the single 
reflector design discussed previously. First, the diffracted is given 
by 


U° IFF = D U 1 ( Q e )i 


-jkp 


/T 


(3.92) 


and 


u 1 * (Q e ) = u* EF 




(3.93) 


which is known from Equation (3.58) and those that follow. Again, one 
obtains the following: 

L ro = p (3.94) 

L rn . PcP (3.95) 

n^ 

p c +p 
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(3.96) 


1 - 1 ♦ 1 

Pi ©i 

c 

E ^ E . 

p i = p c p rl » and (3.97) 


e. 





(3.98) 


Then, the diffraction coefficient is given by 
-j w/4 

D = ~ e - [cot (tt/ 2 + (<fr+i»')/2)F(kL rn a( ♦+♦')) 

2/2iiT 

+ cot(n/2 - ( 4 >+ 4 >‘ )/2)F ( kL ro a ( 4>+4>* ) ) ] (3.99) 


and is known from Equation (3.92). The effect of increasing R c on 

the reduction of the fast varying ripple is shown in Figure 3.18. 

The slow varying ripple may now be reduced by attaching a rolled 
edge to the subreflector (see Figure 3.19). In this case, one finds 


that 


U 


DIFF 

2 


/He gOIFF e ‘ Jkp 2 

/ p c +p2 


where 


DEED 



(3.100) 


(3.101) 


In addition, the range parameters are given by 


L ro = 


DEED 
Pc PI 


DEED . 

Pc + P1 


(3.102) 
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(3.103) 

(3.104) 

(3.105) 

(3.106) 


+ cot ( tt/ 2 - U+n^FfkL^aUn'))] 
and 3” is given by Equation (3.76) and 
S + « = 3"+a 


(3.107) 


(3.108) 


The new diffracted reflected field is then given by Equation (3.100). 

The effect of increasing R c on the reduction of the slow ripple is shown 
in Figure 3.20. 

Another field component that needs to be considered is the 
reflected-reflected-diffracted field (see Figure 3.21). has been 

computed previously and is given by 


58 




59 




DISTANCE FROM CENTER 

(c) R C SUB = 25X. 


Figure 3.20 (Continued). 



Figure 3.21 Reflected -ref lected-di ffracted field 



(3.109) 


,k EDGE 

ir 2 - = U* EF e J P2 


.REF 


Now, the reflected-reflected-diffracted field is given by 


-jkp , DEED 

DIFF RFF e ( Pc P 

U 3 “ U 2 “ 7?" C d (p.*-*',2) + D '“deed - .♦ + * , » 2 )] (3.110) 

Pc + P 

with D(L, $,n=2) as before. From Figures 3.22 and 3.23, one obtains the 
fol lowi ng: 


p E0GE = Fm-Lv + Xs EED - (D E /2) 2 /(4Fm) 
p » C(y-DE/2) 2 +(Lv-x^ EED ) 2 ] 1/2 

f * = it - sin" 1 ( y " D E /2 ) 


(3.111) 

(3.112) 

(3.113) 


a = COS 


i t DEED.2 A 

•1 j ) + 


<pc dee V- 


Fc 


, DEED DEED 
2p i p c 



y = tt/ 2 - e i , 

■ Y + sin" 1 

Now, the diffraction 


and 

PE/2 i . ' 
DEED 

- p i _ 
angles are given 


by 


(3.114) 


e + = 4)+<j>' = e“ + 2 $' , 

. or 

® = 6 + a + 2sin"* 

°E/2 ] 
DEED 

_ p i _ 


(3.115) 
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Figure 3.22 Reflected-reflected-di ff racted field. 
3 . ' e .c 



Figure 3.23 Reflected-reflected-di ff racted field. 
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Note that follows from Equation (3.110). A typical plot Is shown 

in Figure 3.24 that also includes the previously calculated terms. 


As can be seen In the previous figure, the knife edge Is 

undesirable. An elliptical rolled edge Is now attached to the 

subreflector to eliminate the reflected-reflected-diffracted field and 

reduce the diffracted-reflected field at the expense of introducing a 

triple reflected field (see Figure 3.25). This field will be analyzed 
REF 

next. Uj has been calculated previously and 

-Jkp2 


U R 2 EF = U REF e 

where P 2 must be determined. Now the normal of the hyperbolic 
subreflector is given by 


(3.116) 


n = x (a/b )y s y 

/rr(7 s 7b'F 

/l + ((a/b )y s ) 

V 1 + '*s/bF 

and the center of the ellipse is specified by 


(3.117) 


x E x s +Y Eys = x s( D E/ 2 ) x s + (D E /2)y s + Bn 


(3.118) 


(see Figure 3.26). Now the ellipse is paramaterized by 
x = Acosv, and 
y = Bsinv 


(3.119) 

(3.120) 


for 0 < v < 2ir (see Figure 3.27). Tilting the ellipse and shifting its 
center result in Figure 3.28, and one obtains 
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Figure 3.25 Triple reflected field. 
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Figure 3.26 Ellipse addition. 
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Figure 3.27 Ellipse. 



Figure 3.28 Tilted ellipse. 
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x s = xcose - ysin0 + X^, and 

(3.121) 

y s * xsin0 + ycos0 + 

(3.122) 

with _ _ 


e-sm- 1 iE-XsIDt/z) . 
_ 0 _ 

(3.123) 

So the subreflector surface is now completely described. 

Then P 2 and 

IjREF are known. In Figure 3.25, U^^F -j s given by 
2 3 


.REF . .REF r-VT- e J P3 
6 1 / Pc + P3 

(3.124) 

p 3 ■ C(y-y s ) 2 ♦ (lv-x s ) 2 ] 1/2 

(3.125) 

1=2 , and 

p c R c cos 0i 

(3.126) 

R c = l^£A 2 sin 2 v + B 2 cos 2 v] 3 ^ 2 

(3.127) 


where v is given for the ellipse. To determine x s , y s , and 0-j , the 
reflection point on the subreflector is found from the iterative routine 
given in Appendix B. The triple reflected field is then known and a 
plot of this field, the reflected field, the diffracted field from the 
rolled main reflector, and the diffracted reflected field from the 
rolled subreflector is shown in Figure 3.29. 

Two additional fields will now be Included for completeness sake. 
They are the spillover incident field and the reflected field from the 
elliptical rolled edge on the subreflector (see Figure 3.30). Now the 
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with 


1 _ . 1 _ + 2 

Pt f>i R c cosei (3.131) 

and 

R * * [A 2 s1n 2 v + B 2 cos 2 v] 2 ^ 2 • (3.132) 

c AB 

To solve for these variables, the reflection point Is found using the 
same procedure as that In Appendix B. From Figure 3.31 the dot product 
Is given by 

* A 

T »n M r »n 

TTf w 

with 

AAA 

n * v + V 

i * -x - dy , and 

A A 

r * x + fy . 

A 

The normal, n. Is given In Appendix B. Solving for f and d yields 

f = AsinQcosv + Bcoseslnv + YE - t (3.137) 

Acosvcose - Bsinvslne + YE - s 

and 

d s Asinecosv + Bcoseslnv + YE . (3.138) 

Acosvcose - Bsinvsine + XE + Fc-Lv 

Once the reflection point Is known, one finds that 

cose. = Bcosvcose - Aslnvslne + f[Bcosvs1n6 + Aslnvcose] , (3.139) 

[B 2 cos 2 v 4 A 2 s1n 2 v] 1/2 [1 ♦ f 2 ] 1/2 


(3.133) 

(3.134) 

(3.135) 

(3.136) 
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and 


(3.140) 


Pi » [ y r 2 + (Fc-Lv+x r ) 2 ] 1/2 


P r - C(t-y r ) 2 + (Lv-x r ) 2 ] 1/2 . 


(3.141) 


Finally, Figure 3.32 shows the addition of these two additional field 
components as well as those already Included In Figure 3.29. 

As can be seen from this last figure. It appears that the 
Cassegrain system falls miserably when all the major field components 
are included but the last two components examined need not pose any 
concern. First, the feed design Itself will Include a taper so that the 
spillover and single reflected field will be reduced In magnitude. But 
this Is trivial because these fields may be eliminated altogether by 

9 V • 35* 

* 15 * 

« IOOX 

f • eox 

MAJOR AX IS * SX 

9, * 35* MINOR AXIS »2X 



Figure 3.32 Spillover field and reflected field additions. 
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using a pulse radar system. In this system, the length of the pulse 
width determines a window through which the return from the target Is 
examined. Because the pathlengths of the spillover field -and single 
reflected field differ greatly from the desired target return, the pulse 
radar does not see these returns since they are not In that selective 
window. But, the triple reflected field has nearly the same pathlength 
as the desired reflected field so this component can't be overlooked. 

The same Is true with the diffracted fields, but these are under 
relatively good control. So the triple reflected field Is seen to be 
the major obstacle to good system performance. This field will now be 
examined a little more closely. 

Up to this point, the target area or plane of Interest has been 
placed at the focal length. Now It Is advantageous to allow this 
distance to vary for added flexibility (see Figure 3.33). The four 
field components being examined must be modified, and this Is done In 
Appendix C. In Figure 3.34, a representative plot shows the total field 
and the total field less the triple reflected field. The triple 
reflected field at grazing Incidence Is not accurately portrayed in this 
figure. To increase the accuracy, an additional factor is included in 
the triple reflected field expression as follows: 


R 


s,h 



-j£L 3 /12 -jir/4 
e e 


-1 


F(X L ) 


Vu L r 


(3.142) 


But for a magnetic line source 
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*h 


e * J5L3/12 


-j*/4 

e 


-1 

Tipr 


F(X L ) ♦ q*(C L ) 


C L « -2m(Q R )cose^ 

X L - 2kL L cos 2 6 1 

l L « s r s* 
s r +s' 


m( Q) * [k Pg(Q )] 1/3 


(3.143) 

(3.144) 

(3.145) 

(3.146) 

(3.147) 


and Pg(Q) Is the radius of curvature of the surface at the reflection 
point. This Is already known, s r Is P3, and s' Is • so 


L 


L 


P 3 • 


(3.148) 


Note that cos 6^ was calculated before and q* is the Fock Integral which 
Is given by 


q*(«) 


1 J v'(t) e" JgT dT 
77 -• W£ ( t) 


with 


2jv(t) « wj(t) - W£(t) 
and 


w1(t) 

2 


1 

77 


-J i 


s*j 


2ir/3 


Tt - t3/3 

e dt 


(3.149) 


(3.150) 


(3.151) 


and c arbitrarily small and positive. 
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Now inserting this factor yields Figure 3.35 as a representative plot. 

It is apparent that this increased accuracy shows that the triple 
reflected field is a serious problem in the Cassegrain system. 

To begin to reduce this field, first consider a Gregorian 
subreflector system as shown in Figure 3.36. With this offset design 
and the subreflector placed low, the triple reflected field is virtually 
eliminated. But this configuration introduces a doubly reflected field 
that was not present before and must be considered. First, the normal 
reflected field is found (see Figure 3.37). The expression for the 
first reflected field is 


DC _ -J k Pi -jk(p c + Pl) j*/2 

jREF = e yr>c e e 

1 / »i 


(3.152) 


The caustic distance p c may be calculated from geometrical 
considerations. So, one obtains that 


2 2 1/2 

Pc = C ( Lv+x s ) +y s ] 


(3.153) 


2 2 1/2 

Pi = C(Fc+Lv+x s ) +y s ] , and 


(3.154) 


2 1/2 

P 1 = [ ( V y s )2+(Fm+Lv+x s * )3 • (3 ' 155) 

Then U^EF i s gi ven by 
2 

-jkp 2 

e (3.156) 

with 


2 

P2 = DISPLN + Fm+Lv -y m 

Wm 
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(3.157) 




r , - 1 1 — * — i — * ' r 

30. 40. SO. to. 

distance from center 


1 — r 

70. 



Figure 3.35 Total field ( ) and total field less triple reflected 

field ( ). 
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Figure 3.37 Reflected field. 

since the caustic distance in this case is infinity. Finally, y s may be 
expressed in terms of y m by noting that 

TANa = P 2 - DISPLN - Lv = Lv + x s 

^m -ys 


which yields 


/ 


1 + Lv z - 2aLv 
(be ) z 


y s = -Lv + a - a 


(3.158) 


Now consider the doubly reflected field shown In Figure 3.38. Again, 
attach an elliptical rolled edge to the subreflector. The analysis 
proceeds as with the triple reflected field. The angle e needs to be 
examined more closely in this situation. If the ellipse Is attached 
below or at the main axis (Figure 3.39), then 

, . e. • s1n -l[XE- x s (AT T ACH POINT)] 

1 B (3.160) 

but If it is attached above the main axis (Figure 3.40), then 

6 « w - 0j . (3.161) 



Figure 3.38 Doubly reflected field. 
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For a given point In the plane of Interest, two reflection points 
need to be found simultaneously as Illustrated in Figure 3 . 41 . 
Proceeding as In Appendix B, the first dot product is * 



Figure 3.41 Doubly reflected field, 
with 


A w * 

- X _ y Ri y 
n i « . TO 

tt + yin^pr 

4Fm* 


(3.163) 


r l * x + ey, and 


(3.164) 
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Tj ■ x + cy. (3.165) 

Solving for e and c yields that 

y Rl - y (I) (3.166) 

[ y ^” _ - _ (Fm+lv+D I SP L N)J 

4Fm 

and 

c = y Rl * (Asinecosv + Bcosesinv = Y£ . (3.167) 

L y Rl 2 " (Fm+Lv+Acosecosv - BsinQsinv + X^J 
4Fm 

The second dot product is 


A A 


n 2* r 2 a n 2 * r l 



(3.168) 


with 


n 2 a Bcosvcose - AsinvsinQ x + Bcosvsine + AsinvcosQ y (3.169) 
|_B 2 cos 2 v + A 2 sin 2 v] 1/2 LB 2 cos z v + A z sin z vJ i/2 

A A 

F 2 = -X - cy , and (3.170) 

I 2 « -x - ay. (3.171) 


Thus, one finds that 


a _ Acosvsine + Bsinvcose + # (3.172) 

Lv+Fc+Acosvcos 6 - Bsinvsin9 + X£ 

Now the simultaneous solution to the two dot products in Equations 
(3.162) and (3.168) is needed to find the reflection point. This is done 
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Iteratively using bisection methods over the v and 
these values are found within some specified error 
the two reflection points as 


y 


M Intervals until 
R1 

Then v and y give 
R1 



y R1 2 /(4Fm) 


(3.173) 


* Acosvcose - Bsinvslne + Xc, and 
R2 L 

(3.174) 

y R2 * Acosvslne + Bsinvcose + Y£ . 

(3.175) 

Then 


pep -JM -3k pj 

u i F * 6 / p cs e 

1 /pJ ^ PCS+Pl 

(3.176) 

where 


PI * [(Fc+Lv+xr2) 2 + YR2 2 3 1/2 

(3.177) 

p i * C(Fc+Lv+xr2-xri) 2 + yRi - yR2) 2 3^ 2 

(3.178) 

1 1 2 
Pcs s Pi + R c cos6i 

(3.179) 

Rc « 1 (A2sin 2 v + B 2 cos2y)3/2 
W 

(3.180) 

and cos8i may be found from the second dot product. 

Next the second 

reflected field Is 


f p cm e 

2 1 '• Pcm + P2 

(3.181) 

where 


1 * 1 + 2 

Pcm Pirn Rc cos ( 1t "®i) 

(3.182) 
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Pirn " PI + Pcs 


(3.183) 


R = 2Fm 
c 


i + r yR i 2 ) 

{ W ^ 1 


3/2 


COS ( ir— 9-j ) = -COSSi 


and cos9j may be found from the first dot product, 
distance is 


(3.184) 

(3.185) 
Finally the caustic 


P 2 - [(Fm+Lv+DISPLN-y R i 2 /(4Fm) ) 2 + (yRi-y(I)) 2 ] 1/2 (3.186) 

REF 

and Ug is known. Typical plots for the sum of this doubly reflected 
field and the normally reflected field are shown in Figure 3.42. It is 
apparent that attaching the edge higher on the subreflector reduces the 
field ripple but intuitively this will increase the triple reflected 
field effect which was the original problem. In any case, the effect of 
the doubly reflected field with the Gregorian system is unacceptable. 

But this system does provide insight into reducing the triple reflected 
field by using an offset reflector type of design. So, next let us try 
utilizing the offset design with the classical Cassegrain system. 

Initially, this type of design poses problems because edges will be 
attached to the top and bottom of both reflectors resulting in several 
diffracted fields at the junctions. Although the triple reflected field 
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Figure 3.42 Sum of two reflected fields. 
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Figure 3.42 (Continued). 


may be reduced, the additional diffracted fields will degrade overall 
system performance. If these diffracted fields could be reduced and the 
triple reflected field minimized, an acceptable system would be the 
result. Through the blending of the edges considered in the next 
chapter, diffracted fields are reduced sufficiently to allow the offset 
Cassegrain system to work effectively. 
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CHAPTER IV 
THE BLENDED SURFACE 


The last chapter brought to view problems encountered In the 
Cassegrain system. Now a viable solution Is presented to solve these 
problems and obtain satisfactory system performance. 

The subreflector edges are the first that will be blended as shown 
In Figure 4.1. By a blended surface It is meant that the curved 
surfaces are blended between an ellipse and the original surface to 
which these surfaces are attached. As before the hyperbola Is described 
by 


x = a [/I+Iy/bP -1] 


(4.1) 


with 


e 


sin(fc,+$ r )/2 
si n( 4> v -ij> r )/z 


(4.2) 


a 


Fc 

“Ze , and 


(4.3) 


b * a/e z -l 


(4.4) 
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The blended upper surface Is given In general by 


f * f HYPERBOLA P + f ELLlPSE^‘ P ^ 

(4.5) 

P( v ) - 1+cos (irv/vf ) . 0 < P < 1. and 

(4.6) 

v . Vf VJ 


f Vf-Yf , 0 < v < vf . 

(4.7) 

Note that the parameter through which the blended edge 

Is generated Is 

Yp with Y'j , Yf and vf to be specified. Proceeding as 

In the previous 

chapter, the ellipse to be blended Is given parametrically as 

X * Aslnvcose - Bcosvsine + Xell 

(4.8) 

Y = Asinvsine + Bcosvcose + Yell 

(4.9) 

with 


e = si n” 1 f X ELL - X(v=0)^ . 

(4.10) 


Xeh and Yell are calculated as In the previous chapter and Increasing v 
is shown In Figure 4.2. Now Equation (4.5) yields 

X <V ■ *<VHYI>ERBOLA ,, <y + X < Y p >ELLIPS E t ‘"Pf Y p > 3 (4.11) 


Y <y * v p p (y + ''(yELuosE^-xyi <4.121 

with the blending function P given by Equation (4.6) and v given by 
Equation (4.7). 


87 



The bottom blended surface is generated in a similar manner. In 
Figure 4.3, v now varies from a specified v.j to ir/2. The blending 
function is now 


P(v) 




with 


v 



Vi < v < w/2 


(4.13) 

(4.14) 


and Yi and Yf are to be specified. The ellipse to be blended is given 

by 

X = -Acosvcose - Bsinvsine + X^ll (4.15) 

Y = -Acosvsine + Bsinvcose + Ygu_ (4.16) 

with 

- X(v=w/2) 

B J • (4.17) 


9 = ,in- V«ElL 


Again X^y. and Y^l are calculated as in the previous chapter. The 
equations generating this blended surface are then given by 

x ( Y p) = X ( y p )hYPERBOLA^ 1_P ^ Y p^ + VeLLIPSE^V (4.18) 

and 

Y ^V = Y p^ 1_p ( Y p)3 + Y ^ Y p^ELLIPSE P ^ Y p^ * (4.19) 

Of course between the two blended surfaces the normal subreflector is 
described by Equation (4.1). 

Now the centers of the ellipses will be given. The normal to the 
hyperbola is given by Equation (3.117) and evaluated at Y^ for the top 
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edge. Then, one finds that 


X ELL * + YELL Y « X(Yi)X + YjY + Bn (4.20) 

and for the bottom edge with n evaluated at Y f is specified by 

X ELL x + Yell Y * X(Y f )X + Y f y + Bn . (4.21) 

Now the offset subreflector with two blended surfaces Is completely 
specified by Equations (4.18), (4.19), (4.11), (4.12), and (4.1). The 
parameters that can be varied are Yj and Y f for the top and bottom, vj, 
Vf, and the major and minor axes of both ellipses. 

The design of the surfaces may proceed and a few general rules are 
given here to simplify the procedure. It Is convenient to use the same 
values for the axes of both ellipses and fix v-j at -ir/2 and vf at it. 

This allows the knife surfaces to be completely removed from view of the 
source. The excessive hidden surfaces on the subreflector, as well as 
the main reflector later on, are removed by limiting the range over 
which the curve Is generated. This Is In general better than 
respecifying v^ and vf which would lead to surfaces whose shape Is 
changed each time. In general, the larger the surfaces, the flatter the 
field will be. But small surfaces are also desirable to taper this 
field to reduce the effect of the triple reflected field. So a tradeoff 
exists in the size of the surfaces . This size is determined by the 
major axis of the ellipse as well as the length of the blended hyperbola 
section (which is determined by Y^ and Yf). Now the minimum radius of 
curvature of the edge must be kept greater than a quarter wavelength at 
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the lowest frequency of operation. Once the general size of the 
surfaces is decided, the radius of curvature is determined by the minor 
axis dimension of the ellipse. One final constraint is the relationship 
between the major axis length and the length of the hyperbola section. 

By looking at the slope or first derivative of this curve, a gradual 
transition is best obtained for a hyperbola section length that is 
approximately four times the major axis length. Again, this constraint 
is true at Vj equal to -n/2 and Vf equal to it. So when designing the 
surfaces , the size for acceptable performance is the overriding 
consideration. Besides the design of the surfaces, the size of the 
subreflector itself is another variable. It must be large enough to 
illuminate satisfactorily the main reflector yet kept as small as 
possible so that interactions between the two reflectors may be 
minimized. 

The field produced by the subreflector will now be examined. Using 
UTD and physical optics, the field along a parabolic contour is 
calculated. Using a magnetic line source (Figure 4.4), one finds that 
the reflected field is given by 


UREF 


- j k p-j - j k pp 

e /Hie 

/p? / P c +Pr 


(4.22) 


For the region between Y^ arid Yf, U R ^ has been calculated before with 


PC = C(Lv-x)2 + y2] 1/2 


(4.23) 


Pi = [Y2 + (Fc-Lv+x)2] 1/2 


(4.24) 


p r = [( YP-Y) 2 + (Fm-Lv+x-YP 2 /(4Fm)) 2 ] 

90 


1/2 


(4.25) 


c- > 



I SI 

YP = 2Fm [zL + /!_ +1 ] I 4 ' 

c / c 2 

with 


Now the reflected field from the edges is more complicated. Th( 
caustic distance must be recalculated using 


l/p c = 1/ Pi + 2/(R c cos9i) . 

For a given point, p^ is known from Equation (4.24). Now 
R c = | 1 /CUR V 1 


where 


CURV = 


y'x" - x'y" 
[(y') 2+ (x') 2 ] 3/2 


(4 

(4 

(4 


and the derivatives need to be taken with respect to the parameter 
The derivative results for the bottom edge are as follows: 


Y, 


X' = X'hypU-P) + X HY p(-P ' ) + X' ELL P + X ELL P' 

X'HYPW . fy(l+(y/b ) 2 )- l/2 
b 2 

P'(Yp) - - 1 " ) 

2 tt/2-v-i Y f"’i 

X'r|i(i D ) = [Asinvcose - Bcosvsineir y^-vi j 

tLL K Tf-Yi 


(4.31) 

(4.32) 

(4.33) 


(4.34) 
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(4.35) 


X" - X HYP (1-P) 2 X HY p (-P') ♦ Xhyp(-P") 

« I 

+ X ELL P + 2 X ELL P' + X EL L P" 

N 

X HYP( Y P) ■ a/b 2 


[l+(Y/bn 


TTTTT 


P M (YP) - -1 cos 


„ v-vl ' 
W^vT 


YFY7 


* 

Acosvcose + Bslnvsln© 

”ir/2-vi - 



TFTT 


X ELL ( Y p) 

Y ' «(1-P) +Y p (-P') + Y ’ ll P + Y ell P 1 

I 

Y ELL (Yp) * (Acosvsine + Bcosvcose) 


ir/2-Vl“| 

YFYT 


Y" B -2P" ♦ Y p (-P-) ♦ Y« ll P ♦ 2 Y' ll P' ♦ Y ell P” , and 


Y E ii(Yp) s (Acosvsine - Bslnvcose) 


~ir/2-v-j — 

T^YT 


The derivative results for the top edge are as follows: 


*' * X HYP P * W" * X ELlU- p > ♦ * EL l I*') 


P'(Yp) * -1 S!N[*V]( »_ ) 

x ‘ L (Vp) * (Acosvcose + Bs1nvs1ne)( v f ) 

Y f-Yl 


(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 
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X" = X HYP P + 2X hyp P' + X HYp P" + X ELL (1-P) + 2X P1 , (-P 1 ) + X F11 (-P*) 


ELL' 


'ELL' 


* -i cos 

irv 

TT 

7 

Vf 

YT^YT 


X EL l( y p) * (-Asinvcose + Bcosvsine) 


" vf 

YT*T 


Y* • P + Y d P" + Y ell ( 1-P) + Y f , , (-P 1 ) 


ELL' 


■ 

YEL E (Yp) = (Acosvsine - Bsinvcos9) 


" Vf ‘ 


Y" = 2P 1 + Y d P" + Y ell ( 1-P) + 2Yp, , (-P 1 ) + Y f1 , (-P" ), and 


'ELL' 


'ELL' 


(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 

(4.51) 


Y ELL ( Y p) 


(-Asinvcose + Bcosvcose) 



(4.52) 


These derivatives provide the information necessary to compute the 
radius of curvature. The normal is computed next so that cosQ-j may be 
found. For both edges the slope is given by 


m = Y'/X' 


(4.53) 


and the normal has slope -1/m. For both edges the normal is given by 


A A A A A 

n = X x + Y y = -x + (l/m)y . 

[1+1/m 2 ] 1 ^ ( 4 - 54 ) 


Now forming the dot product with the incident vector, one finds that 
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A A * 

Pi , -(Fc-Lv+x)x-Yy 
UI-c- T v ~ +X7*+V*Ji'2 


(4.55) 


which yields 


cos6i « -X n (Fc-Lv+X)-Y n Y 
TT F c - Tv7 X)^y^Ji^ 


(4.56) 


Note that the caustic distance p c Is known. Equating the angle of 

Incidence with the angle of reflection allows the point YP on the 

parabola to be found. So one obtains that 

* * 

p r * -(Fm-Lv+X-Yp2/(4Fm))x+(YP-Y)y 

nTm-l\+X-YP^( 4 hm))^+(YP-Y)^P^ (4.57) 

and 

cosQi = -X n (Fm-Lv+X-(YP) 2 /(4Fm))+Y n (YP-Y) 

C(Fm-Lv+X-Yp2/(4Fm)) 2 +(YP-Y)2]l/2 . (4.58) 

Using an iterative bisection method, YP may be found numerically. Then 
the reflected field from the blended surfaces along a parabolic contour 
is known. Numerically, it Is convenient to pick points on the 
subreflector and ray trace to the corresponding point on the parabolic 
curve. 

The subreflector is now analyzed using the physics optics 
approximation and the results compared with the reflected field given by 
geometrical optics. Again using a magnetic line source, the current on 
the lit side of the subreflector Is given by 

• A m, 

Js « 2n x Hi . (4.59) 
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A A A A A 

t = t x x + tyy = -x-my f or m negative or 

/l +m 2 m positive on lower edge. (4.64) 

Also, one finds that 

Pi = [Y 2 + (Fc-Lv+X) 2 ] 1/2 (4.65) 

and Jt is then known. Now, one obtains the following: 

p = [(Fm-Lv+X-(YPAR)2/(4Fm))2 + (YPAR-Y)2] 1/2 (4.66) 

AAA 2 A A 

p = p x x + V - -(Fm-Lv+X-( YPAR) /(4Fm))x + (YPAR-Y)y 

p (4.67) 
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Figure 4.5 Physical optics analysis. 


SHADOW 


Figure 4.6 Three Integration regions. 
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(4.70) 


Pl « -(Fc-Lv+x)X-Yy . 

Pi 

A 

The normal (n) is known from Equation (4.54) so the remaining dot 
product is known in the integral. The integration is carried out over 
three regions as shown in Figure 4.6. The lit region endpoints 
L 1 and L^) may be found using an interative bisection technique. Once 
these points are known, the integration is carried out suppressing the 
leading constants. Then, the scattered field is given by 

- j2ir(p'j+p) * 

Hz = / (txPy"ty Px) ® ( Pi *n)dYp (4.71) 

Li /pfP 

or numerically 

N ~j 2ir( p-j +p) - - 

Hz * £ (t x Py - ty Px ) 6 ( Pi *n) AYp (4.72) 

i rnr 

where Yp is the parameter defining the curve (or just y in the case of 
the hyperbola section). This is evaluated using the trapezoidal rule 
for numerical integration. 

A typical subreflector is shown in Figure 4.7. The slope or first 
derivative is also present to show the slope transition from the 
reflector to the blended surfaces. The UTD and physical optics plots 
follow. Next, a smaller subreflector is shown with the same size 
surfaces (see Figure 4.8). In this case the subreflector is too small 


♦ v » 45° 

♦r • 5° 

Fm » 81.33X 

Fc - 81.33X 
Major Axis » 2.67X 
Minor Axis * 1.73X 
Hyperbola Section * 10.67X 
Edges attached 
at 4.12X and 6.543X 



(a) Subreflector with blended edges 



Figure 4.7 Subreflector with blended surfaces. 




♦» = 45 ° 

*r * 5° 

Fm = 81.33X 
Fc = 81.33X 
Major Axis = 2.67X 
Minor Axis = 1.73X 
Hyperbola section = 10.67X 
Edges attached 
at 5x abd 5.5X 



(a) Subreflector with blended surfaces. 



(b) First derivative plot. 

Figure 4.8 Smaller subreflector with same size surfaces. 
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in terms of wavelengths. Finally the original size subreflector with 
smaller surfaces is considered (see Figure 4.9). Again, performance 
suffers and the slope transition in this last figure is much more 
abrupt. 


It is also helpful and convenient to look at the far field of the 
subreflector using geometrical optics (GO) and moment method theory. 
For GO, the far field is given by 


-jkpi 

. -J k Ps 


U = e 

/ f»c e 

(4.73) 

/Pi 

/ p c +p s 

* 

ith the variables shown in Figure 4.10. 

Now p s is the far field 

i stance so 

Pc + Ps * Ps 

, and 

(4.74) 

~jk p-j 

-jkp s 


U = e 

/*T e 

(4.75) 

/p? 

/Pi 



But this last factor is a common factor and may be ignored yielding 



(4.76) 


at a specified angle <j> which may be easily found from previous 
calculations. 

For the moment method, points are set up along the curve of the 
subreflector a distance "ds" apart. The distance used was usually about 
0.2 wavelengths at a particular frequency. On the blended surfaces, 

"ds" was approximated as a small line segment. For the hyperbola 
section (Figure 4.11), the uniform segments were found using 
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♦ v « 45° 

♦r * 5° 

Fm = 81.33X 
Fc * 81.33A 
Major Axis * 2. 67 A 
Minor Axis = 0.433 A 
Hyperbola section = 
Edges attached at 
4.12X and 6.543A 



(a) Subreflector with blended edges. 



(b) First derivative plot. 

Figure 4.9 Original subreflector with smaller blended surfaces 
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(4.77) 


dx _ ay/b 2 

*7 ' [l+(y/b) 2 ]i /2 

at the bottom point and 

dy = ds (4.78) 

[l+(dx/dy) 2 ] 1/2 . 

This dy is added to the previous y to obtain the next y from which the 
corresponding x coordinate may be calculated. Once the subreflector is 
represented by the point geometry spaced ds apart, the moment method 
program may be run using this information to obtain far field plots. 
Using the same subreflector geometries as in Figures 4.7 through 4.9, 

GO and moment method far field patterns are shown in Figures 4.12 
through 4.14. The moment method analysis provides accurate results and 
for the relatively small subreflector is easily and quickly implemented. 
These results also show the closeness of including only the reflected 
field in the UTD analysis. When these two procedures yield results that 
differ greatly, then more than the reflected field is being seen, and 
the design must be compensated so that it looks more like a simple 
reflective surface. Once satisfactory performance (flat amplitude over 
area of main reflector to be illuminated) is obtained, the design of the 
main reflector may proceed. 

As with the subreflector, the main reflector will consist of a 
parabolic surface with two blended surfaces attached (Figure 4.15). Of 
course, the parabolic section is given by 

x*y2/(4Fm). (4.79) 
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Figure 
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(b) Moment method plot for geometry of Figure 4.7 (a). 
4.12 Field plots for geometry of Figure 4.7 (a). 
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(a) Geometrical optics plot for geometry of Figure 4.8 (a). 



(b) Moment method plot for geometry of Figure 4.8 (a). 
Figure 4.13 Field plots for geometry of Figure 4.8 (a). 
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(a) Geometrical optics plot for geometry of Figure 4.9 (a). 



(b) Moment method plot for geometry of Figure 4.9 (a). 
Figure 4.14 Field plots for geometry of Figure 4.9 (a). 
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Figure 4.18 Ellipse for bottom edge. 

The blending function used on the top edge is 


cos*v 

P = 1+ 

— 2 


0 < P < 1 


where 


v = Vf Yp_YlT 0 < v < Vf . 

YFT-YIT 

The ellipse in Figure (4.16) is parameterized by 

XE = Asinv, and 
YE = -Bcosv . 

The normal of the parabola is given by 

n = -x+( Y/2Fm) y 
[l+(Y/2Fm) 2 ] 1/2 


(4.80) 


(4.81) 


(4.82) 

(4.83) 


(4.84) 
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so the origin of the blended ellipse is located at 


X ELT 


-B + YIT 2 

cu gn ) 2 ] 1/z , and 

"ZFm 


Y 


ELT = 


B (YIT) 
2Fm 




+ YIT. 


(4.85) 

(4.86) 


the angle of the main axis of the ellipse is now tilted (Figure 4.17) 
and given by 


9 - sin" 1 CYIT 2 /(4FM)-X ELT 3 
B 

so that the elliptical curve is generated by 

X = Asinvcose + Bcosvsine+X^Lj 
and 

Y = AsinvsinQ - Bcosvcose + Y^^j . 

Finally the top blended surface is generated by 

x = + (Asinvcose + Bcosvsine + X ELT )(1-P) 

and 


(4.87) 


(4.88) 

(4.89) 


(4.90) 


Y = YpP + (Asinvsine - Bcosvcose + Yelj)(1-P) (4.91) 

where Yp is a parameter that varies between YIT and YFT. 

For the bottom surface, the ellipse (Figure 4.18) is parameterized 

by 


XE = -Acosv 


(4.92) 
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and 


YE = -Bsinv . 

The origin of the ellipse in this case is at 


X = -B + YFB 2 

ELB [>( YFB ) 7 ] 1/2 
2Fm 


and 


B(YFB) 

Y = 2Fm + YFB 

’elb ii7 - ( - wj2fre 

2Fm 


The tilt angle of the ellipse is 


e = sin" 1 r YFB /( 4Fm )- X ELBl 
B 1 

and the ellipse curve is given by 


X = -Acosvcose + Bsinvsine + X^lb * and 
Y = -AcosvsinQ - Bsinvcose + Y^lb . 

The blending function for this surface is 


P = 



v-v, ■] 

1+COSTT 

I 

ro 

i 

< 


0< P < 1 


with 


= (w/2 - v : l( Y P- YIB )+y. Vi < v < it/ 2 
1 Wfb-ytr ; 1 


"YFB-YIB' 

Then the bottom blended edge is generated by 


(4.93) 


(4.94) 


(4.95) 


(4.96) 


(4.97) 

(4.98) 


(4.99) 


(4.100) 
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( 4 . 101 ) 



( l-P)+( -Acosvcos 0 + Bsinvsine + X^ LB )P 


and 


Y = Yp(l-P) + (-Acosvsine - BsinvcosQ + Y^gjP . 


( 4 . 102 ) 


The entire main reflector surface may now be generated. As in the 
case of the subreflector, larger surfaces lead to flatter field 
performance. This overall size is determined by the major axis length 
plus the parabolic section length blended together. The parabolic 
section length is determined by (YIT,YFT) and (YIB,YFB) for the top and 
bottom surfaces , respectively. The minimum radius of curvature of the 
surface is determined by the minor axis dimension, and this is kept 
greater than a quarter of wavelength at the lowest frequency of 
operation. A final rule of thumb is to make the length of the parabolic 
section about four times longer than the major axis of the ellipse in 
order to make the slope transition from the parabola to the blended 
surface as gradual as possible. As in the subreflector case, v-j is 
fixed at -tt/ 2 , and Vf is fixed at tt. This assures that the surface 
extends well into the back area of the reflector. Since much of this 
blended surface is excessive, some of it may be eliminated by 
controlling the range over which the total reflector is generated. This 
is easier to implement than changing the variables which govern the 
shape of the curve which would result in a different surface design. 

With the geometry given, the field in the plane of interest or 
target area may now be calculated. The main reflector is analyzed alone 
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by placing the source at the virtual feed. The moment method procedure 
and GO are used to analyze the main reflector. First, the moment method 
is used to consider some special cases. 

In Figure 4.19, a typical reflector with sharp edges is shown. The 
moment method plot follows. The ripple may be greatly reduced by 
attaching elliptic rolled surfaces with no blending involved as in 
Figure 4.20. Now using a linear blending function of the form v/vf 
instead of the cosine function described earlier, the field plot in 
Figure 4.21 shows even more improvement. A parabolic blending function 
of the form (v/vf)^ is considered next. The improvement of the field in 
this case is shown in Figure 4.22. Finally, the original cosine 
blending function is used as shown in Figure 4.23. It and the parabolic 
blending show similar results. The cosine blending is chosen since it 
gives the best results and also does a good job of shaping the surfaces 
in the back region. The distance to the target area or plane of 
interest is also a variable, and care must be taken in choosing this 
distance. If the observation plane is too close to the reflector, one 
is limited by the feed position. If the plane is moved too far away 
from the reflector, far field effects become apparent, and the plane 
wave gradually becomes more and more tapered. These effects are shown 
in Figure 4.24 for various observation plane distances where "DISPLN" is 
referenced to zero at the virtual feed. The case for DISPLN * 0 is 
shown already in Figure 4.23b. 
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*v * 45° 

♦r - 5° 

Fm * 610 cm 
Top at 505 cm 
Bottom at 322 cm 


(a) Parabolic reflector section. 



Figure 4.19 Parabolic reflector section and field plot. 
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Top at 505 cm 
Bottom at 322 cm 
Major Axis * 80 cm 
Minor Axis * 30 cm 
Parabolic section 
blended « 320 cm 


ISO. 


(a) Linearly blended edges. 



(b) Moment method plot. Frequency - 4 GHz 
Figure 4.21 Reflector with linearly blended surfaces and field plot 


Figure 4. 



♦ v « 45° 

♦r “ 5° 

Fm * 610 cm 
Top at 505 cm 
Bottom at 322 cm 
Major Axis * 80 cm 
Minor Axis * 30 cm 
Parabolic Section 
blended * 320 cm 


(a) Parabolic blended edges. 
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(a) Cosine blended edges. 
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(b) Moment method plot. Frequency ■ 4 GHz 
Figure 4.23 Reflector with cosine blended surfaces and field plot. 
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The GO reflected field for the main reflector will now be 


calculated (see Figure 4.25). 

In the parabolic region, the reflected 

field is given by 


,jREF . 

/pf 

(4.103) 

with 


Pi = [(Fm-X) 2 +Y 2 ] 1/2 # 

(4.104) 

Pp = Fm + DISPLN-X, and 

(4.105) 

YPLN = Y . 

(4.106) 


For the top and bottom blended surfaces, the reflected field is 


U 


REF _ e 


-jkpi 


/pT 


/ 


~PT~ 

Pc + Pr 


-jkpr 

e 


(4.107) 


with 


1 = 1+ 2 

Pq P'i R c cos ©i (4.108) 


and pi as before. To calcualte R c and cos9i, the partial derivatives 
with respect to Yp must be found. For the bottom surface, one finds the 
following: 


* = *PARA^“^ + *PARfl("P') 


PARA' 


+ x ell p + x ell p* 


(4.109) 
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(a) DISPLN = -305cm. 
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(b) DISPLN = 305cm. 

Figure 4.24 Moment method plots for various distances to the 
observation plane. 
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Figure 4.25 Reflected field from main reflector. 
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( 4 . 110 ) 


X PARA (Yp ) " 





X ELI_( Y P) ■ [Asinvcos 0 + Bcosvsin9](^~^|^) 
X = X PARA^ 1_P ^ + 2X PARA^" p,) + X PARA^” P " ^ 


+ X ELL P + 2X ell P + X ELL P" 


X PARA (Yp) 


1 

‘ZFm 


P" ( Yp) 


= 1 cos 

7 


r v "Vi 
* it/2-vi ^ 


YFB-YIP 


X ELL 


= [Acosvcose - Bsinvsine] 


— / w j %»n 
11/ V I 

YFb-'YIB 


2 


r = (i-p) + Yp(-p-) + y ell p + y ell p* 


Y ELL (Yp) = ( Asinvsin0 " Bcosvcos9)(^£l^_) 


( 4 . 111 ) 


( 4 . 112 ) 


( 4 . 113 ) 

( 4 . 114 ) 

( 4 . 115 ) 

( 4 . 116 ) 

( 4 . 117 ) 

( 4 . 118 ) 
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9 


and 


(4.119) 


Y" = -2P' + Yp(-P" ) + Y^ ll P + 2 Y^ ll P' + Y aL P" 


YeLl(Yp) = [AcosvsinQ + BsinvcosS] 


ir /2-v-j 

YFB-YIB 


For the top surface, one obtains the following: 


(4.120) 


X ' * X PARA P + X PARA P ' + X ELL^ 1_P ^ + X ELL^~ P ^ 


(4.121) 


P' = -I sin [11] f L__) 

2 L Vf n YFT-Yir 


Kn = (Acosvcose - Bsinvsin9)(^^ ) 


(4.122) 

(4.123) 


Y" = X" P + ?X' P + X P" 
PARA PARA X PARA P 


+ x ell( 1“ P) + 2X rn (" p ') + X FM (-P") 


'ELL 


ELL' 


(4.124) 


P"(Yp) = 4 cos [11] 
7 v f 


* “|2 


YFT-YIT 


x EL l( y P) = (Asinvcose - Vcosvsine) 


v f 

YFT-YIT 


Y' = P + YpP' + Y ell ( 1-P) + Y p, , (-P* ) 


'ELL' 


(4.125) 

(4.126) 

(4.127) 


Y' (Yp) = (Acosvsine + Bsinvcos9)f v f ) 
ELL YFT-Yir 


(4.128) 
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9 


and 


( 4 . 129 ) 


Y" = 2P * + YpP" + % L (1-P) + 2 YJ ll (-P') + Y ell (-P") 

Y ELL = ("4sinvsin9 + Bcosvcos 9) (y^y^) 2 
From elementary calculus, the position vector is given by 

A A 

r = xx + yy 

with 

A A 

r' = x'x + y'y» and 

l?'| - C(x') 2 + (y ' ) 2 ] 1/2 . 


( 4 . 130 ) 


( 4 . 131 ) 


( 4 . 132 ) 

( 4 . 133 ) 


Also, the unit vector is 


- r' x'x + y'y 

1 = W\ = C(y') 2 +(y ') 2 ] i/2 


( 4 . 134 ) 


and 


A A 


t » = x [x"-x' (x'x"+y l y")[(x' ) 2 +(y' ) 2 ] -1 ]+y [y"-y' (x'x M +y 'y")[(x‘ ) 2 +(y* ) 2 ]~*] 

( 4 . 135 ) 


[(x') 2 +(y') 2 ] 1/ " 


Ci 1 u 
l 111 01 I I JT j 



The radius of curvature is given by 



( 4 . 136 ) 


( 4 . 137 ) 
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and the outward facing normal by 


n =11 


(4.138) 


The dot product with the incident vector 


Pi = (Fm-x)x-yy 
[(Fm-x) 2 +y 2 ] 1/2 


(4.139) 


rvtrkw 4* U U ~ -C J i - 1 -I J - 

majr men uc iui mcu jneiumy 


COS 0-j = pj •n 


(4.140) 


and the caustic distance p c is known. 

The reflected field is solved completely when YPLN and p r are 
known. First, one obtains that 


P r * C(Fm + DISPLN - x ) 2 + (YPLN - y ) 2 ] 1/2 


(4.141) 


For the bottom edge (see Figure 4.26), one finds that 


,-l 


* TAN ( (Fm-x)/y) 


and 


-1 


02 = cos (cose^) . 


(4.142) 

(4.143) 


The normal is always outward facing so if > ir/2 then subtract a 2 from 
it to get desired 02 * Then, one obtains that 


TAN( ai +2o2) = ( Fm + DISPN " x )/ (y — YPLN) 


(4.144) 


or 
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YPLN = -(Fm + DISPN - x)/TAN(o 1 + 2^) + y 

Similarly for the top surface, one obtains that 

TAN(<* 1 + 202 - tt/ 2) * (YPLN - y)/(Fm + DISPN - x) 


(4.145) 


(4.146) 


YPLN = (Fm + DISPN - x) TAN(a^ + 2 - ir/2) + y . (4.147) 

The reflected field for the surfaces is then given by Equation (4.107). 

In Figure 4.27, a typical main reflector is shown. The slope or 
first derivative is also present to show the slope transition from the 
reflector to the blended surfaces. The GO and moment method plots 
follow. Next, smaller surfaces are attached (Figure 4.28) and the 
degradation in the field is apparent. The GO plot is not as greatly 
affected since it only contains the reflected field components. Larger 
surfaces will result in flatter field performance, but a tradeoff must 
be made between edge size and field performance. Once the field is 
acceptable, the surfaces should not be increased further in order to 
keep the triple reflected field to a minimum. 

The total system may now be put together. The subreflector and 
main reflector have been analyzed separately, and there is little 
further design that can be done at this stage. The offset between the 
reflectors should be as large as possible, and the blended surfaces have 
been made as small as possible to reduce the triple reflected field 
without sacrificing flat field performance. 
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(c) Geometrical optics plot. Frequency = 4 GHz 



(d) Moment method plot. Frequency = 4 GHZ 
Figure 4.27 (Continued). 
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(a) Main reflector with blended edges. 



(b) First derivative plot. 

Figure 4.28 Parabolic reflector with smaller blended surfaces. 
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IfO. 209. 256. 909. 362. 915. 966. 521. 579. 627. 690. 

DISTANCE FROM AXIS (CM) 


(d) Moment method plot. Frequency = 4 GHz 
Figure 4.28 (Continued). 
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Moment method and GO analysis are used for the entire system and 
the GO reflected field for the entire system will be examined first. 

In Figure 4.29, the reflected field from the subreflector that does not 
interact with the blended surfaces is given by 


D rr -J‘ k (pi+Prl+Pr?) 

U REF = / Pci e 

v Pi (Pcl + Prl) 

(4.148) 

with 


p d = [(Lv-x $ ) 2 + y $ 2 ] 1/2 

(4.149) 

and 


Pi = [(Fc-Lv+x s ) 2 + y s 2 ] 1/2 . 

(4.150) 

Now y m and y s are related by 


y m = 2F m [-l/c + /17c 7 + 1] 

(4.151) 

where 


c = ^s . 

Lv-x s 

(4.152) 

Also, one finds that 


v = w 2 / f nc~\ 

"m ^m / v-»« »•/ 

/ A 1 CO \ 
\ -r • J- -J'J J 

p r 2 = Fm + DISPN - x m 

(4.154) 

y ob = ym » and 

(4.155) 

2 2 1/2 
Prl » C(Fm-Lv+x s -x m ) + (y m -y s ) ] . 

(4.156) 
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Points on the subreflector are then chosen and rays traced. The 
hyperbola section of the subreflector reflects rays that also extend 
outside the parabolic section of the main reflector. The reflection 
point on these surfaces must be found differently. 

The reflected fields for the top and bottom surfaces in Figure 4.30 
are given by 


y REF 


-jkpi 

e 

/pT 


/ 


PcT~ 

Pcl+Prl 


-jkprl 

e 


/~P cT~ 
'f P c2 +, V2 


-jkpr2 

e 


(4.157) 


or 


U REF 


. -jk(pi+Prl+Pr2) 

/ PclPc2 e 

* Pi (Pcl + Prl)(Pc2 + Pr2) 


and 


(4.158) 


2 2 1/2 

Pi = [(Fc-Lv+Xg) + y s 3 


(4.159) 


where x s and y s are chosen on the subreflector. There are two 
possibilities for p c i. If (X S ,Y S ) is on the hyperbola section then 

Pci = [(Lv-x s ) + y s ] (4.160) 


Otherwise, p c i for the surface s is given by Equation (4.28) and those 
that follow it. The reflection point on the main reflector must now be 
found. For (X S ,Y S ) on the hyperbola, one obtains 


n i = V + v = 


A 

-X + 


1 + 


(a/b2)y S - 
[l+( Ys/b) 2 ] 1/2 x 
(a/b 2 )y s ~ 
[l+(ys/b) 2 ] 1/2 I 


2 1/2 


(4.161) 
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In the subreflector section, cose^ was calculated such that 


cose. = nj*p rl 


and 


P r i = -(Fffl-L v +y s -ym)x+(y m -y s )y 
Prl 


where 


r i 


= r ( Pm- 1 v/4-v * - 
''b 


'Mil/ 


( w . 

\«/m 


/ 


. 2 -. 1/2 


(4.162) 

(4.163) 


(4. 164) 


for a given parameter Yp for the main reflector. The reflection point 
( x m>ym) tlien found using a bisection routine to solve for Yp. For 
the blended surfaces, the normal was calculated in the subreflector 
section, and the reflection point is found in a similar manner. The 
normal for the main reflector was also previously calculated, and the 
incident vector is 

A A 

p i = “ p rl ‘ (4.165) 

So, one finds that 

cose. = p..n 2 (4.166) 


and P r 2 and YOB may be calculated as in the main reflector section. 
Finally, the caustic distance is given by 


!_ « 1 + 2 


Pc 2 Pi RcCosQj 

(4.167) 

with 


p i = p rl + p cl 

(4.168) 
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and R c was calculated previously. Then i$ known through selecting 
points on the subreflector and implementing the above ray tracing 
procedure. 

In the moment method procedure, the entire system is represented by 
a geometry of points spaced a variable distance apart. This distance is 
usually fixed at two-tenths wavelengths for best results. Now the 
moment method is applied to the entire system (Figure 4.31) and the plot 
in Figure 4.32 results. The magnetic line source also reflects off the 
main reflector into the target area, and this field is subtracted 
yielding Figure 4.33. Next the spillover incident field and the direct 
reflected field off the subreflector edge are subracted out yielding 
Figure 4.34. The triple reflected field ripple is apparent as well as a 
slow varying ripple caused by successive reflections from the 
subreflector edge, the main reflector edge, the subreflector hyperbolic 
area, and the main reflector parabolic area. Of course in an actual 
system, this field component is not seen by using a pulsed radar system. 
Finally, by selectively zeroing the correct elements in the impedance 
matrix, the interaction between the main reflector and subreflector is 
eliminated resulting in Figure 4.35 where the slow and fast ripple have 
been eliminated. The corresponding GO plot is shown in Figure 4.36. 

GO may also be used to calculate the triple reflected field level 
since all the basic analysis has already been completed (see Figure 
4.37). First °2 is recalculated and (x s ,y s ) found on the subreflector 
using a bisection method. Then P r 2 is known. The normal is computed at 
(x s ,y s ) and another bisection method to find YOB. This yields P r 3 and 
p c which may also be calculated at this point. 


140 



4V = 45° 

4>r = 5 

Fm = 610cm 
F c = 610cm 

Main Reflector Subreflector 


Major Axis = 80cm 

Minor Axis * 30cm 

Parabolic blended section = 320cm 

Top edge at 505cm 

Bottom edge at 322cm 


Major Axis = 20cm 

Minor Axis = 13cm 

Hyperbolic section blended = 80cm 

Top edge at 49.4cm 

Bottom edge at 30.9cm 


Unnecessary edge extensions in back regions have been eliminated. 


Figure 4.31 Entire Cassegrain system. 
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Figure 4.33 Entire system less main reflector contribution. 
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Figure 4.34 Entire system less main reflector, subreflector and 
spillover contribution. 
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Figure 4.35 Desired reflected field. 
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Then 
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(4.169) 


A plot of this field is shown in Figure 4.38. The level of this field 
is plotted relative to the reflected field magnitude in the center of 
the target zone. The triple reflected field is actually higher than 
that which would be given by the maximum level in Figure 4.38. This 
could be remedied by using a more accurate value for the reflection 
coefficient when computing the triple reflected field. 

At this point, little else can be done about the triple reflected 
field without sacrificing the desired flat field performance. One 
alternative is to place microwave absorber in the vicinity of the 
subreflector to block the ray path of the triple reflected field from 
the subreflector to the target area. This improvement could be easily 
implemented. In actuality, the triple reflected field will be weaker in 
a three-dimensional system in that additional spread factors will reduce 
the triple reflected field even further. 

Finally, the source horn is placed at the "origin" for convenience 

WPl 1 nntimal norfnrmanro Paco in rhanninn hnnnc anri horn 

” ' ' ’ ' ~ 1 * •• 

orientation and the close proximity of radar components make this an 
ideal location. A typical broadband horn gives the desired flat field 
over angles much greater than those considered here. If possible, a 
field taper is desired over the blended surfaces which will improve 
system performance but the major purpose of the horn is flat, broadband 
performance over the widest possible angles. 
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Figure 4.38 60 triple reflected field. 
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CHAPTER V 


CASSEGRAIN SYSTEM DESIGN PROCEDURE 


The Cassegrain system is a viable alternative for providing a 
uniform plane wave in a compact range measurement system. A procedure 
will be described briefly for designing a Cassegrain system given a 
target area dimension as a constraint. No further constraints will be 
given but in actuality each individual application will have its own 
unique restrictions. 

The system described in Chapter IV will be used as an example; that 
is to say the target area is to have a dimension of six feet. The main 
reflector is considered first. It is convenient to make the vertical 
dimension the same six feet. Increasing this dimension will flatten the 
field when just considering the main reflector, but this will result in 
a larger overall size leading to an increased triple reflected field. 
Therefore, the main reflector is kept as small as possible to keep the 
triple reflected field to a minimum. The overall focal length (Fm) is 
chosen to be approximately three times the target dimension. Care must 
be taken to make sure the target plane is not in the far field of the 
main reflector at the minimum frequency of operation. For this 
particular case, the focal length was chosen to be twenty feet. 
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The edges of the main reflector must be considered next. Only the 
field of the main reflector (source at virtual focal point) is examined. 
The three blended surface variables are the major axis length, minor 
axis length and blended parabolic section length. The blended parabolic 
section is fixed at four times the major axis length for the most 
effective transition from the parabola to the blended surface . The 
major axis length should be increased until the field over the target 
area is flat. The natural taper of the field will also be present. 

Some rounding of the field at the junctions will be tolerable since the 
subreflector field is usually tapered and will compensate for this 
problem resulting in a smooth field at the junctions. Finally the minor 
axis length is usually one-third to two-thirds the major axis length. 

The actual length is set such that the minimum radius of curvature of 
the edge is greater than a quarter wavelength at the lowest frequency of 
operation. The rounding of the field at the junctions for this case had 
a 0.1 dB to 0.2 dB variation as is seen in Figure 4.27d of the previous 
chapter. 

The subreflector is designed next. The goal is to generate a 
uniform field illuminating the parabolic section of the main reflector 
with a subreflector whose overall dimensions are kept as small as 
possible to reduce field interactions between the two reflectors. One 
parameter that determines the size is $ r . For this example a value of 
five degrees was chosen. The blended surfaces are added next 
following the same procedure as with the main reflector. The far field 
pattern is then examined and a half dB variation over the area to be 
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illuminated is tolerable. In Figure 4.12b the subreflector field plot 
is shown and the desired illumination area is from 135° to 150° for this 
particular example. If the field variation is not acceptable, the size 
of the blended surfaces may be increased. If a flatter field is not 
obtained, the size of the hyperbolic reflector section must be 
increased, and the addition of the blended edges repeated until a 
satisfactory pattern is obtained. 

The total system may then be put together, and the resulting 
pattern examined. The feed is placed beneath the main reflector at the 
"origin". Although the feed could be placed somewhere between the 
reflectors, this location is most operationally convenient and gives 
good source performance since the beamwidth is smaller and field flatter 
at this increased distance. The system is now completely specified, and 
little further design may be done except for varying the offset angle 
between the two reflectors. Increasing this offset angle without 
altering the rest of the system should reduce the triple reflected 
field. The interaction between the blended surfaces of the reflectors 
now results in weaker reflected fields at each edge yielding an overall 
reduced triple reflected field. This also changes the desired field 
pattern somewhat since the positioning of the hyperbolic subreflector 
and parabolic main reflector has changed. This change usually does not 
affect the desired pattern greatly but the triple reflected field is 
reduced. Therefore, the ideal situation is to have a maximum offset 
angle to reduce the triple reflected field making sure that the desired 
field over the target area remains acceptable. 
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For this example, consider three cases with offset angles of 35°, 
45°, and 55°. A system with = 35° is shown in Figure 5.1. The 
minimum triple reflected field from Figure 5.2 is about 1.6 dB. The 
desired reflected field with interactions between reflectors eliminated 
is shown in Figure 5.3. Within the six foot target area there is about 
0.5 dB rolloff at the extremes. There is about 0.2 dB variation over 
4.8 feet of the area. The next case in Figure 5.4 has <t> y = 45° and was 
the system considered in the previous chapter. At this angle the triple 
reflected field is reduced to about 0.9 dB (see Figure 5.5). The 
desired reflected field has a 0.4 dB rolloff within the six foot area 
and 0.2 dB variation over 5.4 feet (see Figure 5.6). The final case has 
<j>v = 55° as shown in Figure 5.7. The triple reflected field has now 
been reduced to a level of 0.6 dB (see Figure 5.8). This change of 0.3 
dB is not as great as the reduction between the angles of 35° and 45° 
which implies diminishing returns with increasing angles. The desired 
reflected field has less than 0.4 dB rolloff within the six foot target 
area. And 0.2 dB variation over 5.4 feet as shown in Figure 5.9. 
Therefore, a maximum offset angle is desired but this must be tailored 
to each individual case. Note that the polarization performance will 
deteriorate with increasing tilt angle for a full three-dimensional 
system. 
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CASSEGRAIN SYSTEM (CM) 


*v = 35° 
V - 5° 

Fm = 610 cm 
Fc = 610 cm 


Main Reflector 


Subreflector 


Major Axis = 80 cm 

Minor axis = 30 cm 

Blended parabolic section = 320 cm 

Top edge at 384 cm 

Bottom edge at 201 cm 


Major Axis = 20 cm 
Minor Axis * 13 cm 
Blended hyperbolic Section 
Top edge at 47.4 ccm 
Bottom edge at 24.5 cm 


Unnecessary edge extensions in back regions have been eliminated 


Figure 5.1 Cassegrain system with ^ = 35°. 
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Figure 5.3 Moment method plot of desired reflected field. 
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C95SEGBRIN SYSTEM (CM) 


<frv « 45° 

*r * 5 ° 

Fm = 610 cm 
Fc = 610 cm 


Main Reflector 


Major Axis = 80 cm 
Minor axis = 30 cm 
Blended parabolic section 
Top edge at 505 cm 
Bottom edge at 322 cm 


= 320 cm 


Subreflector 

Major Axis = 20 cm 
Minor Axis = 13 cm 
Blended hyperbolic Section 
Top edge at 49.1 ccm 
Bottom edge at 30.9 cm 


Unnecessary edge extensions in back regions have been eliminated 
Figure 5.4 Cassegrain system with <j> v = 45 0 . 




Figure 5.5 Moment method plot showing triple reflected field. 



Figure 5.6 Moment method plot of desired reflected field. 
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CASSEGRAIN SYSTEM (CM) 


549 . 610 . 


<tV = 55° 
4>r = 5° 

Fm = 610 cm 
Fc = 610 cm 


Main Reflector 


Subreflector 


Major Axis = 80 cm 

Minor axis = 30 cm 

Blended parabolic section = 320 cm 

Top edge at 635 cm 

Bottom edge at 452 cm 


Major Axis = 20 cm 
Minor Axis = 13 cm 
Blended hyperbolic Section 
Top edge at 50.29 cm 
Bottom edge at 35.35 cm 


Unnecessary edge extensions in back regions have been eliminated 
Figure 5.7 Cassegrain system with 4 > v = 55°. 
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8 Moment method plot showing triple reflected field. 
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CHAPTER VI 


CONCLUSIONS 


In conclusion, the Cassegrain system will provide a uniform plane 
wave but the triple reflected field ripple must be reduced through 
judicious design and absorber blocking. The blended edges provide 
superior performance over large target areas for a given size reflector 
system. Design and analysis is best implemented by studying the main 
reflector and subreflector separately. Then the system may be put 
together to complete the analysis. 

Different techniques were used to study the system. The moment 
method provides accurate results but is limited by structure size and 
computational speed. Though difficult to implement in some situations, 
UTD provides results which compare favorably with moment method even 
when just examining the reflected field. UTD is also fast and usable on 
large structures. Both these techniques are best utilized by examining 
each reflector spearately before combining and checking the total 
system. 

With the two reflector system, the elimination of undesired field 
components becomes the prime consideration. Those field components 
which have pathlengths that differ greatly from the desired reflected 
field pathlength may be eliminated through time gating with the use of a 
pulsed radar system. The triple reflected field component does have a 
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similar pathlength to the desired reflected field and must be reduced 
through careful design of the two reflector system. If 
three-dimensional effects were also considered, the apparent level of 
this field would be reduced further. Absorber blocking around the 
subreflector would also reduce the triple reflected field at the expense 
of introducing diffracted components from the absorber. The diffracted 
fields from all blended surface junctions have been virtually eliminated 
through the blending process. 

This report is not a complete study of this topic, but it serves to 
illustrate the potential benefits and problem areas associated with 
subreflector compact range systems. Further work is needed for 
three-dimensional structures in order to give additional insight into 
the working of the system as well as give improved accuracy in the area 
of several field components. The blending procedure implemented is new 
and various blending functions as well as different blending processes 
provide additional possibilities for system enhancement. The diffracted 
field at the blended surface junction, though not considered here, is 
another area of potential analysis. Other Cassegrain reflector system 
configurations are possible including the isolation of one reflector 
from the other in separate areas to reduce the interaction between the 
two. Finally, the actual physical construction and implementation of a 
Cassegrain reflector system would provide the final verification of the 
analysis and design. 
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APPENDIX A 


REFLECTION POINT ON MAIN REFLECTOR 


This appendix describes how the reflection point on the main 
reflector is found. Knowing that the angle of incidence equals the 
angle of reflection, a point, (X 0 ,Y 0 ) is initially chosen on the 
parabola (see Figure A.l). Now the normal is given by 


n a x _ 




y 

2 yr rz 


A new coordinate system is formed at "0" (see Figure A. 2). 
the radius of curvature and given by 


3/2 


R = 2f(l + y ) 
TF 7 


(A.l) 


R is 


(A.2) 


so "0" is located at 


X 0R ’ X o + 2f d + 
and 


(A. 3) 



Working in this new coordinate system, move A<f> and assume R remains the 
same for this small change. Now 

COS0j = COS 0p 
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or 

A A A A 

n • i = n • r 

and 

A A A A 

T = (x'x + y'y) - xRcosa<|> - yRsinA<|> , 

A A 

r = x(x" - Rcos a$) + y(y" - Rsi n A<|>) , 

and 

A A A 

n = -(xcosA<t> + ysinA<|>). 

Now as A<{> 0, cos A<|> + 1 and sin A<t> A4>. So 

T - x(x'-R) + y(y'-RA<|>), 

r - x(x"-R) + y(y"-RA<|>), 


(A. 5) 


A A A 

n - -(x+yA4>). 

So Equation (A. 5) simplifies (ignoring second order terms, (A<j>)2) to 

A , . (x l '-R) 2 [(x l -R) 2 + (y') 2 ]-(x t -R) 2 [(x ll -R) 2 + (y H ) 2 ] . 

2(x'-R)y 1 [(x"-R) 2 +(y") 2 ]-2y"R(x'-R) 2 +2y 'R(x"-R) 2 -2(x"-R)y" • 

[(x' -R) 2 +(y ' ) 2 ] 

(A. 6) 

Now moving from this nearby point to a point actually on the reflector, 
the error is computed as follows: 

A A A A 

e = n NEW* 1 NEW -n NEW' r NEW * ( A * 7 ) 
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The reflection point is found when the error is below some prescribed 
value. Otherwise, the procedure is repeated until the minimum error is 
obtained. 

For convenience, the point actually on the reflector is shown in 
Figure A. 3 and given by 

Yp — Yqp ^ Rsin(A<f>+a) 

and 

Xp = Yp 2 /(4f) 

where 

a = TAN-l(Y 0 /2f) . (A. 8) 



.11 m . 

t (x , y ) 
PLANE 


Figure A. 3 Actual point on reflector. 
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Then 


A A 

x-lty 

(y *OLD " y POLD^ y » 

A 

^ y "oLD ‘ y POLD^ y * 

The subscript "OLD" refers to the original coordinate system. The error 
is then given by Equation (A. 7) as 


NFW = Zf 

(1 + Yp 2 /4f 2 ) 1/-2 


T NEW = 'oLD “ x P0LD^ x + 


and 


'NEW = ( x "0LD " x P0LD^ x + 


G = 


iTTiWj 


rrr 


/(X ‘oLD ‘ X P0LD* Z + ^ Y 0LD " Y P0LD^ 


Y 

fX‘ -X -fY 1 _ Y ^ POLD) 
1 OLD *P0LD 1 OLD ’POLD'-gf — 


A x “old " X P0LD^ 2 + ^ Y "0LD ‘ Y P0LD^ 


(X" 


OLD “ X P0LD 


(Y" 


OLD 


Y P0LD 


jYPOLD) 


(A.9) 


Finally, the points (x',y*) and (x",y") need to be transformed to 
the new coordinate system as follows (see Figure A. 4): 


X ' NEW = (x 0R " x 'oLD^ COSot ‘ (y 0R ” y 'oLD^ sino 


(A. 10) 
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(a) (x',y') coordinate transformation. 


A 



(b) (x",y") coordinate transformation. 
Figure A. 4 Coordinate system transformation. 
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y, NEW * " C(x 0R ‘ x '0LD )sina + {y 0R * (A. 11) 

and 

X "nEW = ^"oLD ” y OR) sina “ ^ x "oLD " x 0R^ COSa (A. 12) 

y "NEM * ( y “0L0- y 0R )coSl1 + (X "0LD ' *OR )si " a * (A.13) 

» 

Another method of proceeding is to take the cross product 

A A A A 

i x n = n x r 

and form 

A A A A 

n • i _ n • r 

A A A A 

i x n n x r 

or 

AAA A A A A A 

(n • i )(n x r) = (n • r)(i x n) 

or 

A A A A 

(n • T)(n x r) = (n • F) (T x n) . (A. 14) 

A 

Using the small argument forms of T, F, and n. Equation (A. 14) yields 


A<(> 


y' (x"-R) + y"(x'-R) 

( v " -R \ + y" fv'.fti' _ ~7U T u VC 

\ •' ■’/ - \” •*/ -V J 


(A. 15) 


when ignoring second order terms. Equation (A. 15) may be used as an 
alternative to Equation (A. 6). 
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APPENDIX B 


REFLECTION POINT ON ELLIPTICAL EDGE OF SUBREFLECTOR 


This appendix describes how the reflection point on the elliptical 
rolled edge of the subreflector is found. In Figure B.l point (s,t) is 
known and the normal to the surface of the ellipse is given by 


^ A A A <A 

n = x o x + y 0 y = Bcosvcose - Asinvsine x + Bcosvsine + Asinvcose y . 

[B 2 cos 2 v+A 2 sin 2 v] 1/2 [B 2 cos 2 v+A 2 sin 2 v] 1/2 


The dot product is then formed 


(B.l) 


T • n _ F • n 

with 

A A 

T = x + fy , 
and 

* 

r = -x . 

Now at point (s,t), 

t = fs+g or g = t-fs 
so 

y = f(x-s) + t 
or 

f = y^t (B.5) 

x-s 

where for the ellipse 


(B.2) 

(B.3) 

(B.4) 


y = AsinGcosv + Bcosesinv + 
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(B.6) 



and 


x = AcosvcosS - Bsinvsine + Xg . (B.7) 

to find v which gives the reflection point, an initial guess is made. 
From Equation (B.5) f is calculated and then the dot product in Equation 
(B.2) calculated. The difference is formed 

A A 

ERROR = 1 * n - r * n (B .8) 

|T| |r| 

and is this behaves nicely, v may be found by successively bisecting the 
v interval until the desired error tolerance is obtained. Once v and 
the reflection point are known, remaining parameters may easily be 
calculated which finally yield the desired field. 



Figure B.l Reflection point on elliptical rolled edge. 
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APPENDIX C 


MODIFICATIONS FOR VARIABLE DISTANCE TO PLANE 

When making the plane of interest a variable distance (see Figure 
C.l) , the following modifications must be done* For the reflected 
field, p r 2 becomes p r 2 + DISPLN. For the diffracted field from the main 
reflector rolled edge junction, p must be recalculated. First 

9, ■ cos- 1 r w J , (C.l) 

ofl + p c 

0 = tan' 1 (DISPLN) (C.2) 

£ y ( i ) 

and 

p = C(p rl + p c )2 + Y(I)2+ displn2 - 2 (Prl + P^)( Y (n 2 + DISPLN 2 ) 172 

1/2 

cos(e 1 + e 2 )] . (c.3) 

Now 

B" = (C.4) 


or 


6 - 


cos-1 


-y ( I ) 2 -DISPLN 2 + p 2 + 

(Prl * p c > 2 

2p(p rl + p c 



(C.4) 
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Next 


and 

and 

So 


*+ 4 >' 

1/2 sin 

ir/2 - 0 


3 + = 2 <(>* + 




3 + = it - sin“*( g" 1 / 2 ) + 3 ". 

pE,+pt 

K rl v c 

For the diffracted field from the subreflector edge junction, 

x" = Fm + DISPLN. 

OLD 


Finally for the triple reflected field, s goes to Lv + DISPLN 
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